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INTRODUCTION 


A.  General  and  Historical 

Vulcanized  natural  rubber  is  characterized  by  an  ability  to  undergo 
large  reversible  deformations.  In  the  19th  Century,  Couth  ^0,  Kelvin 
[2]  and  Joule  [3]  studied  the  peculiar  thermodynamic  effects  which  acom- 
panied  deformation,  e.  g.  ,  the  reversible  heat  evolution  on  stretching 
and  the  positive  temperature  coefficient  of  the  force  necessary  to  maintain 
a  given  elongation.  These  properties,  as  we  know  now  are  not  unique  to 
natural  rubber,  but  are  common  to  all  crossllnked  elastomers. 

In  the  classical  theory  of  elasticity,  applicable  to  small  strains,  the 
behavior  of  an  isotropic  elastic  body  under  any  type  of  deformation  can 
be  described  completely  if  two  material  constants  are  known,  one  related 
to  the  response  to  shear,  the  other  to  volume  change.  However,  since  this 
theory  is  based  upon  the  limiting  condition  of  Hooke's  law,  it  is  not  applicable 
to  rubber-like  materials  capable  of  large  elastic  deformations.  A  continuum 
theory,  developed  in  particular  by  R.ivlin  t4l  rseks  to  establish  in  general 
terms  the  symmetry  relation  which  governs  the  behavior  of  elastic  materials 
under  different  types  of  deformation.  By  comparing  the  resulting  general 
relation  with  the  experimental  stress-strain  behavior,  one  hopes  to  deduce 
the  equivalent  to  Hooke's  law,  in  terms  of  a  finite  number  of  material  constants. 

One  might  hope  also  to  deduce  such  a  law  from  our  molecular  picture 
of  rubbers  as  an  entangled  mass  of  randomly  coiled  molecules,  i.  e. 
deriving  energetics  from  a  statistical -mechanical  point  of  view.  The  compli¬ 
cations  in  this  approach  are  that  the  picture  of  rubbers  as  above  is  really  that 
of  a  very  complex  viscoelastic  liquid  and  that  the  approximations  used  tc 
derive  an  equation  of  state  on  a  molecular  basis  must  be  severe. 

It  is,  of  course,  also  possible  to  proceed  empirically  to  find  an  applicable 
form  of  Hooke’s  law  by  starting  from  phenomenological  observations  and 
casting  them  into  functional -analytical  form  by  observing  the  necessary 
conditions  of  thermodynamics  and  mechanics.  This  leads  to  workable  equations, 
without,  however,  permitting  a  mechanistic  understanding. 


Ill  -  z 


B,  Treatment  of  Network  Elasticity  in  the  Gaussian  Approximation 

By  1932  Meyer,  Susieh  and  Valko  [5]  ,  Mark  and  Guth  [61  and  Kuhn  [?1 
had  suggested  that  the  elastic  retractive  force  of  rubbers  ari  ses  from  the 
decrease  of  entropy  associated  with  the  smaller  number  of  conformations 
available  when  this  network  of  chains  is  deformed  by  stretching.  The  re¬ 
tractive  force  of  a  single  Gaussian  chain  can  be  very  easily  derived  • 
hence  most  of  our  present  theoretical  knowledge  concerning  rubber-like 
elasticity  stems  from  treatments  invoking  the  Gaussian  chain  approximation 
and  extending  it  to  the  entire  network.  Applying  this  approximation  simplifies 
obtaining  the  changes  in  conformational  entropy  of  the  network  chains,  but 
leaves  a  number  of  difficu'ties  associat'd  with  the  "liquid-like"  nature  oc 
bulk  elastomers  unresolved,  such  as  the  volume  requirements  of  the  molecules 
themselves,  the  volume  dilation  due  to  the  isotropic  component  of  the  applied 
tension,  the  free  energy  change  accompanying  this  volume  change,  and  other 
possible  free  energy  changes  which  might  arise  from  a  change  of  shape  at 
constant  volume.  Many  of  the  early  treatments  assumed  simply  that  all 
these  effects  could  be  neglected  in  comparison  to  the  contribution  from  the 
change  in  conformational  entropy.  Recent  experimental  work,  however  has 
prompted  a  reexamination  and  modification  of  earlier  theories  [8,  9l  . 

The  first  step  in  the  molecular  theory  of  rubberlike  elasticity  is  the 
derivation  of  the  statistical  properties  of  a  single  chain.  It  is  necessary  to 
Know  the  free  energy  of  the  molecule  quantitatively  as  a  function  of  its  end- 

to-end  distance  r  and  the  relation  between  a  given  distance  and  the  force 
necessary  to  maintain  it  .  Because  of  the  Brownian  motion  of  the  chain,  the 
equilibrium  distance  r  must  be  considered  to  be  an  average  value  and  the 
same  must  be  true  for  the  end-to-end  distance  during  the  period  of  application 
of  a  force  f  .  Assuming  true  randomness  of  the  sequential  arrangement  of 
segments,  the  probability  of  a  given  end-to-end  vector  from  r  to  r  /  dr  in  the 
absence  of  an  external  force,  W(r),  can  be  approximated  b/  the  Gaussian 
distribution  function: 


W(r)  dr  :  (9/^T)^  exp  (  -  ^  r^ 


4 ffr2  dr  (1) 
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where  the  parameter  hai  the  value  if  2r0  (rQ  li  the  mean-square 
end-to-end  distance  of  the  polymer  molecule  unperturbed  by  volume  effects). 
The  error  incurred  by  using  the  Gaus*an  function  In  Equation  (1)  has  been 
found  to  be  negligible  provided  r  is  less  than  one-half  of  the  fully  extended 
chain  length  f  10^  .  The  energetics  are  introduced  by  means  of  statistical 
thermodynamics,  relating  W(r;  and  the  Helmholtz  free  energy  of  a  set  of 
chains.  By  differentiation  with  respect  to  r  one  arrives  at  f  ,  the  force 
needed  to  maintain  a  specified  vectorial  length  r: 

f  ;  3  k  T  r  /  (2) 

The  steps  of  this  derivation  include  considering  that  the  entropy  of  a  chain 
will  decrease  upon  stretching  as  the  chain  assumes  various  conformations 
which  are  characterized  by  average  greater  values  of  r  (the  entropy  must  be 
a  maximum  at  the  most  average  probable  value  of  r)  .  From  the  function  of 
the  Helmholtz  free  energy,  A,  in  the  form  of  the  conformational  integral,  the 
retractive  force  is  derived  as  the  tangent  with  respect  to  L  as  given  by 
Equation  (2). 

The  second  step  in  the  theory  consists  in  calculating  the  elastic  free 
energy.  Ae|  ,  of  the  network  as  a  function  of  the  macroscopic  parameters 
which  characterize  the  deformation.  The  network  can  be  considered  as  an 
ensemble  of  chains,  each  "chain"  now  being  that  portion  of  an  original 
macromolecule  which  extend  from  one  crosslink  to  the  next.  All  these  chains 
together  will  determine  the  retractive  force  of  the  crosslinked  sample,  and 
the  behavior  of  this  system  is  assumed  to  be  the  sum  of  the  contributions  of 
v  chains.  Each  chain,  characterized  by  r^~  in  the  undeformed  state  and 
by  r  in  the  deformed  state,  contributes  according  to  Equation  (2).  Since  the 
material  is  supposed  to  be  initially  isotropic,  one  assumes  that  the  components 
of  the  vectorial  length  of  each  chain  are  changed  by  the  external  deformation 
in  the  same  ratio  as  the  corresponding  dimensions  of  the  sample.  This  is 
known  as  an  "affine"  deformation.  It  is  then  possible  to  write  the  expression 
for  the  retractive  force  in  the  case  of  a  simple  elongation  as: 
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V  k  T 


U-*'X2>  (5) 


where  v:  number  of  chains  between  crosslinks,  k:  BoUsmann's  constant, 

Li  :  original  length  of  sample,  1  :  L/Lj  :  elongation  of  sample,  L: 
length  of  deformed  sample,  and  ~r^  /  *r|'  is  the  so-called  front  factor,  in 
which*r^i  is  the  mean-square  end-to-end  distance  in  the  undeformed  isotropic 
stats  of  a  network  for  all  the  chains  having  the  same  contour  length,  and 
7^  is  the  corresponding  mean-square  end-to-end  distance  for  the  undeformed 
chains  without  crosslinks.  Hence,  the  parameter  (jr2  /  r“  )  '  *  represents 

the  geometrical  mean  of  the  linear  deformation  existing  in  the  undeformed 
network  relative  to  a  state  where  the  end-to-end  distance  for  the  set  of  chains 
coincides  with  that  for  the  chains  unrestricted  by  the  intermolecular  links. 

For  the  purpose  of  this  study  it  will  be  helpful  to  reiterate  some  of  the 
more  important  assumptions  and  approximations  used  in  the  various  derivations 

1.  The  v  chains  of  the  network  are  long  enough  (the  degree  of  crosslinking 
not  excessive)  and  volume  effects  negligible  in  order  to  assure  the 
validity  of  the  Gausaian  distribution. 

2.  The  initial  macromolecules  are  of  very  high  molecular  weight  in 
order  to  assure  a  negligible  percentage  of  free  ends. 

3.  The  contributions  by  the  chains  to  the  total  retractive  force  are 
additive. 

4.  The  deformation  of  the  chains  must  be  af'ine  to  the  deformation  of  the 
sample,  which  takes  place  at  constant  volume. 

As  stated  these  approximations  are  rather  restrictive  and  tend  to  over¬ 
simplify  many  features. 

Recently,  the  consequences  of  one  of  the  major  assumptions  underlying 
the  simple  Gaussian  treatment  -  namely  that  the  chains  are  volumeless  and 
do  not  hinder  each  other  -  have  been  reappraised  [111-  Taking  the  volume 
of  the  chains  into  account  effects  the  conformational  entropy  of  the  network. 
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In  the  unstrained  state,  the  presence  of  inter  segmental  obstruction  serves 
only  to  modify  the  effective  step  length  of  the  Brownian  motion  but  does  not 
affect  the  form  of  the  distribution  function.  However,  when  the  chains  are 
stretched,  the  medium  is  no  longer  Isotropic  and  the  end-to-end  distribution 
for  a  particular  chain  becomes  affected  by  its  chain  volume.  Jackson,  Shan 
and  McQuarrie  ['ll] take  into  account  the  change  of  conformaftonal  entropy  by 
intermolecular  obstruction  by  modifying  the  end-to-end  distribution  function 
in  what  they  call  a  self-consistent  way.  If  the  usual  "kinetic-theory"  expression 
can  be  represented  by: 

f  prop  (X  -  X  “2  )  (4) 

then  the  equation  of  state  for  rubber  elasticity  of  these  authors  can  be 
represented  by: 


f  prop.  (\-\  "2)  [1  -  _ ( 2 y  4  X-1)  -A  .  .  .  ]  (5) 

3M(1  -  fs) 

where  f  andX  have  the  same  meanings  as  in  the  previous  expression,  but 
where  f8  is  the  fraction  of  space  occupied  by  polymer  molecules,  and  M  is 
the  number  of  equivalent  statistical  chain  segments.  Plotting  f  /  (X  -X_Z) 
as  a  function  of  (ZXZ  A  1/X)  should  give  a  linear  relationship  with  a  constant 
slope.  This  was  confirmed  by  using  Roe  and  Krigbaum's  data  for  Viton  A  [18] 

Taking  a  usual  value  of  M  -  100,  fB  ,  the  fractional  occupied  space  turned 
out  to  be  0.  94  which  appears  to  be  reasonable.  The  authors  conclude  that 
"One  can  at  least  assert  that  the  volume -exclusion  of  network  chains  accounts 
for  part  of  the  term  of  the  Mooney-Rivlin  equation.  Other  effects,  such 
as  chain  entanglements,  network  connectivity,  and  nonequilibrium  states  may 
possibly  be  the  other  contributing  factors"  . 


C.  Thermodynamic*  of  Rubber  Elasticity 


Considering  a  reversible  process,  application  of  the  first  and  second 
laws  of  thermodynamics  gives 

d  E  =  T  dS  dW  (6) 

where  £  is  the  internal  energy,  S  the  entropy,  T  the  absolute  temperatur 
and  W  is  the  work  done  by  the  surroundings  on  the  system.  If  P  Is  the 
external  pressure,  V  the  volume  and  f  the  force  of  extension,  then 

d  W  —  fdL  -  p  u  V  (?) 

Writing  the  first  law 

dE  =  T  dS  4-  fdL  -  p  dV  (8) 

and  introducing  Gibb's  free  energy,  F  :  H  -  TS  with  H  -  E  -A  p  V 


dF  -  dE  4-  P  dV  4- 

V  dF  ‘  TdS 

-  SdT 

(9) 

dF  =  -  SdT  -A  V  dp 

J.  fdL 

(10) 

or 

(8F/SL)  =  f 

(ID 

P.T 

From  F  “  H  -  TS, 

(rt  F  /  SL)  =  (*H  /  *L)  -  T(SS  /  3L)  (12) 

p.T  p,  T  p,  T 


and  equating  Equations  (11)  and  (12), 

f  =  (>H/*L)  -  T(dS/3L)  (13) 

Pi  T  p.  T 


m  - 

Equating  further  the  second  partial  derivatives  of  K  with  respect  to  T  and  L, 

(*f/*T)  :  -  (*S/*L)  (14) 

p.  L  p, L 

wc  obtain  the  thermodynamic  equation  of  state  of  the  rubber, 

f  =  (*H/*L)  4-  Tf&F/aT)  (15) 

p,  T  p,  L, 


These  quantities  are  all  at  constant  pressure  and  temperature  which 
corresponds  to  most  experimental  conditions.  However,  the  quantities  of 
real  interest  are  the  partial  differential  expressions  at  constant  volume  and 
temperature.  At  constant  pressure  the  positive  change  in  volume  of  the 
rubber,  referred  to  as  dilation  or  dilatatlon{increased  molecular  separations) 
gives  rise  to  a  change  in  enthalpy.  This  change  in  addition  to  that  of  the 
entropy  presents  a  condition  precluded  b  y  constant  volume. 

To  obtain  the  equation  of  state  at  constant  volume  and  temperature, 
again  the  Helmholtz  free  energy  is  applied 

A  =  E  -  TS  (16) 


A  derivation  analogous  to  the  above,  and  remembering  that 


TP  f/3T)  =  -  T(aS/*L)  (17) 

V,  L  V,  L 


leads  then  to: 


f  =  (*E/dL)  4-  T(*f/eVT)  (18) 

V,  T  V, L 

The  problem  can  be  reformulated  by  dividing  the  total  retractive  force, 
f  into  the  two  components 


i  1  £.  -A 


(19) 
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where 


f  :  retractive  force  due  to  change  in  internal  energy, 

e 

f  :  retractive  force  due  to  change  in  entropy 

Past  authors  have  frequently  neglected  the  energy  contribution  to  the 
retractive  force  Tl2,  13] at  moderate  strains.  But  later  work  has  shown 
that  fe/f  can  be  an  appreciable  fraction.  For  an  ideal  rubber,  for  which 
the  internal  energy  does  not  depend  upon  the  elongation  r/rj  of  the  chains, 
the  force  f  is  proportional  to  the  absolute  temperature  T.  However,  if 
the  internal  energy  of  the  chains  depends  upon  r,  then  not  only  r^/r^ 

— 9  O 

but  even  r  will  depend  upon  the  temperature  and  the  force  f  will  no  longer 
be  strictly  proportional  to  T. 

Having  stated  that 


f  =  fe  -A  T(*/3T^  L 


we  can  write 


or 


=  f  -  T<*Wy,L 

!  =  -  T2  ,  SR /TV 

e  '  _  ) 

*  T  V,  L 


(20) 


(21) 


(22) 


from  which  jne  can  derive  that 


fe/f  -  -  T  ^ln[f/T]  j 

&T  V,  L 


(23) 
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Previously  we  found  that 


f 


VkT 


(x -xz)  D 


or 


f/T  - 


function 


-1 


(3) 


and  therefore,  with  equation  (14),  we  find: 


£e/f  =  -  T  (  ainCf/T]  )  -  T(*lnr*o  ) 


ST 


V,  L 


9  T 


V,  L 


(24) 


An  important  advantage  of  the  preceding  derivation  is  that  the  theory 
of  rubber  elasticity  offers  the  possibility  of  interpreting  the  energetic 
component  fe  in  terms  of  the  molecular  properties  of  the  chain. 

An  important  correlation  can  now  be  established  between  the  elastic 
behavior  of  a  crosslinked  amorphous  polymer  and  the  properties  of  the 
same  polymer  in  dilute  solution.  At  the  theta-temperature,  the  intrinsic 
viscosity  of  a  polymer  having  a  molecular  weight  of  M  can  be  expressed 
as  [271 


[  W]  *  ^  (r2  /  M)3/2  M1/2  (25) 

B  0 

where**  Is  a  constant  independent  of  temperature.  Consequently 

2/3  (^In  [t?  /  d  T)  =  d  In  r2  /  dt  (26) 

0  o 

and  from  Equations  (24)  and  (26) 

f./f  =  2/3  ’  TdlnTfjl  /  9T 

e  6 


(27) 
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This  equation  establishes  the  interesting  possibility  to  determine  the 
ratio  of  the  energetic  to  the  entroplc  contribution  of  rubber  elasticity 
from  the  variation  of  the  intrinsic  viscosity  for  various  solutions  of 
varying  theta  temperatures  or,  vice  versa,  the  latter  from  force-temp¬ 
erature  measurements. 

The  calculation  of  f#/f  via  stress-strain  measurements  necessitates 
measurements  at  constant  volume  which  are  difficult  to  perform.  However 
differentiation  and  rearrangement  of  Equation  (24),  gives  £  14,  15]  : 

(a  in  [f/Tl  /  aT)  -  -  (d  in  r2  /  dT)  -  3  *  /  (X3  -  1)  (28) 

Pi  L  ° 


and  hence 


Vf 


SlnCf/T] 

T  ( - —  )  -  3  ffT  /  (X  -  1) 


(29) 


where  ®is  the  linear  expansion  coefficient  of  the  polymer. 

The  sign  of  f  will  depend  upon  the  sign  of  d  in  r2  /  dT,  where  the 

o 

dependence  of  r  upon  temperature  is  determined  by  the  potential  associated 
with  internal  rotation,  ^or  a  long  chain  made  up  of  bonds  with  uncorreiated 
rotations  which  are  hindered  according  to  one  or  another  type  potentials, 

1  mmy 

the  sign  of  d  in  r  /  dT  is  easily  predicted.  If  the  potential  of  internal 

bond  rotation  for  a  given  polymer  is  such  that  the  deepest  minimum  is  near 

the  transposition,  then,  r"2  will  decrease  by  increasing  temperature,  f 

o  e 

will  be  negative  and  the  stretching  process  is  favored  by  a  decrease  of 
internal  energy  or  by  lowering  T.  If  however  the  rotational  minimum  is 
far  from  the  transconformation,  T2  may  Increase  with  increasing  temp¬ 
erature,  since  many  of  the  chains  which  were  in  conforamtions  become 
longer  by  increasing  the  vibrational  energy  and  going  to  transconformations. 
Then  f  will  be  positive  and  the  stretching  process  will  be  opposed  by  an 
increase  in  the  internal  energy.  Only  in  cases  where  all  the  possible  minima 

~"'7 

were  isoenergetic  will  r  be  temperature  independent  and  f  be  equal  to  zero. 
“  o  e 

A  decrease  of  entropy  will  always  oppose  the  stretching  process. 
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There  are  several  different  methods  of  calculating  fe/f  for  a  given 
specimen.  If  the  measurements  can  be  carried  out  at  constant  volume 
then,  by  application  of  Equation  (17),  it  would  give  us  an  immediate  measure 
of  fQ.  However  direct  measurements  remain  difficult,  only  a  few  data  exist 
r?l  .  If  the  change  in  volume  were  known  as  a  function  of  elongation,  one 
could  calculate  f  via  Equation  (17)*  by  applying  : 


(*f/*T) 

V,  L 


(*f/*T) 


P.  L 


(*V/*L) 

p,  I 


(*V/*T)p>l 

fav/*)T,  L 


(30) 


Unfortunately  though  the  cubic  thermal  expansion  coefficiently /ST)  ) 

and  the  compressibility  (  (?>V/8p)-,  L  )  are  readily  measured  as 
long  as  the  rubbers  are  isotropic^  the  change  of  volume  upon  elongation 
is  very  small  (of  the  order  of  10~  cc/cc)  and  quite  difficult  to  measure. 
Holt  and  McPherson  !"l6"l,  and  more  recently  Hewitt  and  Anthony  [l7l 
reported  the  only  available  figures  in  the  field. 


The  calculation  based  on  Equation  (29)  supposes  that  the  rubber  network 
obeys  Gaussian  statistics.  Rearranging  Equation  (29)  gives  us 


,  .  3 

Vf  =  1  -  T(-ii)  -  3  «  T/  (X  -  I)  (31) 

T  T  P,L 

The  relationship  is  exact  but  is  based  upon  the  Gaussian  network  theory.  In 
contrast,  the  presence  of  nonvanishing  values  of  in  the  literature  indicates 
considerable  deviation  of  the  real  network  from  the  simple  Gaussian  theory, 
and  therefore  Equation  (31)  cannot  be  strictly  justified.  Accepting,  however, 
that  C  arises  from  mostly  failure  to  attain  equilibrium  conditions,  then  only 

M 

C j  in  Mooney's  equation  represents  an  equilibrium  property,  with  the  rubber 
obeying  the  Gaussian  theory 


f  :  2Ci  U-X'2> 


(32) 


from  which 


fe/f  =  1  -  T/Cj  •  (aCj/a T)  4-  rvT 


(33) 
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Also  in  this  case  the  values  of  f  /f  calculated  by  application  of  (31)  or 

(33)  should  be  equal.  Checking  with  literature  data,  shows  that  this  is  not 

so  f8^(Fig.  1  ).  Similarity,  forcing  the  data  to  comply  to  Equation  (32)  by 

setting  C  -  0,  should  show  iji  to  be  independent  of  X  .  which  is  not 

2 \ 

the  case.  Even  more  disturbing  is  the  fact  that  the  approximations  employed 
in  the  Gaussian  treatment  should  become  more  exact  as  X  approaches  1.  yet 
it  is  just  in  that  region  that  fe/f  varies  the  most. 

Looking  closely  at  fe/f,  r  1 8l  ,  (Fig.  2).  as  a  function  of  X  >  one  can 
see  that  the  contribution  of  internal  energy  is  real,  is  c**en  quite  large,  and 
can  be  positive  or  negative.  The  strong  dependence  of  fQ/f  on  X  runs  counter 
to  the  prediction  from  the  Gaussian  theory  while  the  drop  in  fe/f  at  large  , 
observed  for  natural  rubber,  could  be  associated  with  stress  induced 
crystallization  by  Smith,  Greene  and  Clferri  [19^  ,  the  rapid  variation  of 
f  /f  at  very  low  X  i®  surprising  since,  from  the  nature  of  the  approximations 
employed,  the  Gaussian  theory  would  be  expected  then  to  be  obeyed  more 
exactly.  A  decrease  in  fe/f  ,  which  was  observed  with  cis-1,  4  polybutadiene 
[20]  ,  and  similar  trends  can  be  detected  in  other  systems,  appears 
to  he  a  general  phenomenon.  There  is  as  yet  no  satisfactory  explanation  from 
a  molecular  point  of  view.  Most  probably,  this  phenomenon  is  associated 
the  value  of  C2.  Further,  the  sign  of  C ^  seems  to  be  always  positive  which 
appears  to  correlate  with  the  fact  that  fe/f  always  decreases  withX  ,  ir¬ 
respective  of  whether  fe/f  is  positive  or  negative.  Splitting  the  term 
into  its  own  entropic  and  energetic  components  shows  that  about  50%  of  the 
magnitude  of  C2  arises  from  energy  effects  in  natural  rubber  [8]  .  Finally, 

(.  contributes  proportionately  more  to  the  total  stress  at  low\  ,  showing 
a  *.ln  that  the  deviation  from  the  Gaussian  theory  is  more  pronounced. 
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D.  Phenomenological  approach  to  Rubber  Elasticity 

For  this  approach  one  must  first  understand  the  nature  of  the  concept 
of  the  stored  energy  function  W,  which  is  defined  as  a  function  of  the  strain 
parameter  \  ,  and  of  the  thermodynamic  parameters  specifying  the  environment, 
such  as  the  temperature  T  or  the  pressure  p.  The  function  W  as  a  form  of 

free  energy,  represents  the  total  amount  of  work  stored  elastically  in  the 
body  in  the  process  of  deforming  it  into  the  state  of  strain  represented  by 
\ .  W,  is  furthermore  a  function  of  state  depending  only  on  the  Initial  and 
final  free  energy  levels.  Theories  based  on  W  are,  therefore,  concerned 
with  equilibria  analogous  to  classical  thermodynamics,  and  complications  must 
arise  from  non-equilibrium  processes  when  comparing  experimental  data 
with  the  theory. 

Since  the  stored  energy  does  not  depend  on  pure  rotation  of  a  body,  W  is 
a  function  of  a  symmetric  matrix,  X  '  ,  representing  a  pure  distortion  of  the 
shape.  A  suitable  rotation  of  the  coordinate  axes  reduces  this  matrix  to  a 
diagonal  form,  whose  three  elements  are  the  principal  extension  ratios 
(\{)  =  (\i  A  2  ,  X3  )  as  defined  earlier.  If  the  body,  as  assumed,  is 

isotropic  in  the  undeformed  state,  we  may  choose  as  the  independent  variables 
for  W  some  functions  of  the  elements  of  which  are  invariant  with  a  rotation 
of  the  coordinate  axes.  There  are  a  large  number  of  conceivable  strain 
invariants  from  which  sets  of  three  must  be  selected  which  are  mutually  in¬ 
dependent.  The  simplest  choice  of  such  a  set  is  given  by 


X  2  £.  X  2  +  X  2 

1  a  3 


(34) 


A  *  z  4-  »■  ,*  1  4-  x  2zx  J2  (35) 
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which  ara  even,  symmetric.  functions  of  X  ^  .  X  ,  and  X^  .  Physically, 
ths  third  Invariant  1^  is  rslatsd  to  ths  volume  deformation  by 


I-  :  (V/V  )  2 
3  o 


(30 


where  V  and  V  are  the  volumes  of  the  bodv  before  and  after  deformation 

Q 

respectively  . 

The  stored  energy  W  as  a  function  of  the  strain  invariants  can  be  ex¬ 
pressed  in  terms  of  a  power  series  t4]  : 


W  - 


00 

£ 

U.  k  r  0 


Cijk  <*1 


3)1  (I, 


3)j  (I, 


-  D 


(38) 


where  the  coefficients  C^k  are  functions  of  T  and  p  only.  An  approximation 

to  W  can  be  obtained  by  retaining  only  a  finite  number  of  low  order  terms. 

If  only  the  first  term  is  retained,  we  have 


WI  :  C100  (I1  ’  3>  (39) 

which  is  an  expression  equivalent  to  that  obtained  from  the  statistical 
theory  of  rubbers  when  approximating  the  rubber  chains  by  Gaussian  coils. 
If  the  second  term  of  Equation  (38)  is  also  retained, 


W2  :  C100  ^2  "  3)  -A  Cqjq  (I2  -  3)  (40) 


which  is  identical  with  the  Mooney-Rivlin  function  for  incompressible  rubbers 
f4]. 


In  the  case  of  simple  elongation,  and  assuming  incompressibility: 


X  \  :  A  «  L/L0 


(41) 
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where  and  L  are  the  lengths  of  the  material  In  the  original  and  stretched 
state,  respectively.  Deformation  in  the  lateral  direction  is  then  given  by 


x  i  :  x  3 


(V/X  vo) 


1/2 


(42) 


The  principal  stresses,  (force  acting  in  the  direction  of  the  principal 
strain  X  ^  ,  and  measured  per  unit  area  of  the  deformed  body),  are  obtained 
from  W  by: 


t 


i 


<Xi/!3 


>Tw/v0] 

( - -  ) 


p 


(43) 


Replacing  values  for  X  ^  andX  3  •  in  Equation  (43)  ,  solving  for  p  and 
substituting  for  the  value  of  p  in  the  expression  for  t,  obtained  from  Equation 
(43)  ,  we  have  that 


f  -  2(-^  •  -y  )  (X  -  1/X2)  (44) 

f  is  now  the  force  of  elongation  per  unit  cross-sectional  area  of  undeformed 
material  [4l  . 

When  the  stored  energy  function  is  reduced  to  the  first  term  Equation 
(39)  becomes: 

«  :  2C100(I^T73->  »«> 

ijoVo 

comparison  of  Equation  (45)  with  the  results  of  the  statistical  theory  of 
Gaussian  networks  then  indicates  that  Rivlin's  constant  Cjqq  is  given  by 

C100  =-fvkT(^)2/3  W 

2  V 

* 

where  v  is  the  number  of  network  chains  in  the  specimen,  k  is  Boltzmann's 
constant,  T  is  the  absolute  temperature,  and  V*  is  the  reference  volume 
Independent  of  pressure  |"18^ 
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(Kg. /cm?) 


Fig,  3  ,  Ref  fl8]  :  Biaxial  stretching,  shear  and  simple 
elongation  (dashed  line  represents  the  prediction  of 


Gaussian  theory) 
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Treloar  Tl2^  determined  the  stress-strain  relationship  for  vulcanised 
natural  rubber  in  aimple  elongation,  uniform  biaxial  elongation  and  for 
atrip  biaxial  etretch.  Qualitatively  he  could  recognise  two  typee  of  deviation 
from  the  prediction  of  the  Gaussian  theory  which  requires  that  ?  W/^I_  -  0, 

M  * 

(  Fig.  3).  firstly,  at  large  strains,  the  stresses  tend  to  Increase  much  more 
rapidly  than  expected  from  the  theory.  Secondly,  even  at  moderate  strains 

the  curves  derived  from  different  types  of  strain  do  not  coincide  with  each 
other  ,  which  indicates  that  *  W/rt  1^  -  l  0.  At  large  strain,  the  be- 

havlor  could  be  expressed  only  by  a  large  number  of  terms  with  arbitrary 
coefficients  for  Equation  (39). 

On  the  other  hand,  the  upturn  in  stress  at  large  strains  has  been  suc¬ 
cessfully  explained  in  terms  of  the  finite  extensibility  of  the  network  chains 
r  1 2*1  .  In  this  region  of  strain,  where  an  appreciable  proportion  of  the  chains 
becomes  highly  exiended,  the  Gaussian  statistical  treatment  is  no  longer  valid. 
We  quoted  earlier  that  departures  from  Gaussian  statistics  first  become 
important  when  the  vector  length  of  the  chain  is  between  one -third  to  one- 
half  of  its  fully  extended  length.  ^211.  In  the  case  of  a  network  of  chains,  in 
which  the  chain  length  is  related  inversely  to  the  degree  of  cross-linking,  this 
implies  that  the  Gaussian  theory  becomes  increasingly  inadequate  as  the  degree 
of  cross-linking  is  increased.  At  the  limit  a  ,  "non-Gaussian"  theory  becomes 
essential  even  under  the  smallest  strains. 
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E.  The  Problem  of  C 

Behavior  at  moderate  strain  strongly  suggests  that  the  Mooney-Rivlin 
function,  Equation  (44)  ,  is  a  better  approximation  than  the  single-term 
Gaussian  function,  Equation  (45)  .  When  one  is  interested  in  the  stress - 
strain  relationship  in  simple  elongation  only  it  is  widely  recognized  that 
the  Mooney-Rivlin  function  is  a  very  good  representation  of  the  data.  This 
holds  not  only  for  natural  rubber,  but  also  for  a  wide  variety  of  rubberlike 
materials  such  as  crosslinked  silicones  and  polyethylene,  polyisobutylene, 
and  Vitonfluoro  elastomers. 


No  generally  acceptable  molecular  interpretation  of  the  C  term  has 

2 

emerged  yet.  Gaussian  statistics  predict  a  zero  value  of  ,  but  most 
experimental  work  shows  the  existence  of  a  C  term  as  discussed  before. 

As  will  be  shown  later  ,  extremdy  precise  measurements  are  needed  in  order 
to  determine  the  value  of  C2  at  small  strains. 


Values  of  C  may  be  of  the  same  order  of  magnitude  as  C  .  In  a  few 
2  1 
cases  is  greater  than  Cj  ,  but  there  are  many  cases  reported  where  C 

is  very  small.  There  is  no  clear  cut  rule  by  which  one  could  predict  the 

magnitude  of  from  knowledge  of  the  structure  of  the  material.  The  situation 

is  different  for  the  behavior  of  which  increases  regularly  with  increasing 

crosslink  density  and  with  increasing  temperature,  as  follows  from  the 

Gaussian  theory.  Adherence  to  the  latter  could,  as  stated  earlier,  require 

that  the  C2  term  goes  to  zero  at  low  strains.  Van  der  Hoff's  [22^  very 

recent  data  tends  to  show  that  this  might  in  fact  be  the  case. 


Experiments  on  swelling  of  vulcanized  rubbers  show  that  deviations  from 
the  Gaussian  theory  become  less  noticeable  with  dilation  [23,  24]  i.  e.  with 

reduced  rubber -rubber  interactions.  In  the  work  of  Gumbrell,  Mullins,  and 
Rivlin  r24'  ,  the  effect  of  various  swelling  agents  on  the  stress-strain  re¬ 
lationship  of  natural  rubber,  butadiene -styrene,  and  butadiene -acrylonitrile 
copolymers,  showed  that  in  all  cases  that  the  Mooney-Rivlin  2  term  equation 
was  well  obeyed  and  that,  if  the  force  per  unit  cross-sectional  area  of  the 
unswollen  specimen  is  expressed  as  : 

2  1/3 

f  :  (X-X  *)  (2C*  ■/-  2C'  A.  )/v  _  (47) 

u  12^ 

(where  V2  is  the  volume  fraction  of  the  rubber.  )  C'j  ,  the  equivalent  of  C 
for  the  swollen  network,  1 
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becomes  independent  of  V£,  while  the  value  of  C1^,  the  equivalent  of 

C  for  the  swollen  network,  decreaees  Unearily  with  decreasing  v  ,  and  be* 

2  2 
comes  ze*o  at  v  ^  _  o,  25  irrespective  of  the  nature  of  the  solvent.  C'^  can 

therefore  be  represented  by 

C'  ;  C  (4v  -  1)  (48) 

2  i  2 

3 

where  is  the  value  of  the  constant  for  the  unswollen  rubber.  Other  workers 
have  essentially  verified  these  results  [25]  with  crosslinked  networks  of 
silicone  rubber,  butyl  rubber,  and  polymethyl  methacrylate. 

Another  important  observation  was  made  by  Ciferri  and  Flory,  i.  e.  that 

the  magnitude  of  can  be  influenced  by  the  time  scale  of  the  measurements 
[25]  .  Their  samples  exhibited  considerable  stress  relaxation  when  held  at 
constant  length  and  the  apparent  values  of  Cg  decreased  by  10  to  40  percent 
when  the  time  of  relaxation  between  successive  elongations  was  increased 
from  15  mins  to  10  hours.  A  similar  observation  of  the  effect  of  time  was 
made  by  Kraus  and  Moczvgemba  [26]. 

Thus  the  values  of  C may,  at  least  partly,  be  affected  by  failure  to  achieve 
equilibrium  conditions  during  stress -strain  measurements  or  when  insufficient 
time  is  allowed  for  completion  of  stress  relaxation.  But  as  decay  of  stress  in 
stretched  elastomers  is  often  approximately  linear  with  the  logarithm  of  time 
[24,  27,  281  ,  it  is  difficult  to  know  when  the  criterion  of  equilibrium  has  been 
adequately  fulfilled.  This  led  Ciferri  and  Flory  [25]  to  suggest  that  under  ideal 
equilibrium  conditions  the  values  of  C  might  be  expected  to  become  negligibly 

_  M 

small.  Blatz  and  Ko  [29J  reported  that  for  polyurethane  rubbers  which 
exhibited  no  detectable  relaxation,  the  value  of  C£  was  zero. 

Other  facts,  however,  lead  one  to  suspect  that  the  equilibrium  value  of 
C2  may  not  vanish.  Swelling  the  samples  with  a  solvent  and  then  deswelling 
them  improves  their  approach  to  equilibrium  [29]  .  When  this  technique  was 
used  by  Ciferri  and  Flory  f25*1  in  a  butyl  rubber  vulcanizate,  the  value  of  2C 
was  found  to  decrease  to  0.  5  kg/cm4  from  the  value  of  1.  20  kg/cm  obtained 
without  the  swelling  and  deswelling  cycle,  but  did  not  become  zero.  Roe  and 
Krigbaum  [8,  30"1  studied  vulcanizates  of  natural  rubber  and  Viton  floroetastomer 
by  allowing  relaxation  for  at  least  24  and  48  hr,  and  yet  the  values  of  2Cg 
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obtained  were  0.  88  and  2.  37  kg/cm2  ,  respectively.  These  vatues  of  2C^ 
are  certainly  not  negligible 

Contradictory  results  are  also  reported  on  the  relationship  between  the 

magnitude  of  C2  and  the  crosslink  density.  Gumbrell,  Mullins  and  Rivlin  [24] 

prepared  a  series  of  sulfur  vulcanisates  of  natural  rubber  with  2C  •  2.  0 
2  2  5 

kg/cm  ,  although  their  values  of  2Cj  ranged  from  2.  0  to  6.  2  kg/cm  Kraus 

and  Mocsvgemba  C  26]  show  that  when  polybutadiene  potymers  of  different 

primary  molecular  weights  were  crosstinked  with  controlled  amounts  of  sulfur 

so  as  to  give  approximately  the  same  value  of  Cj  ,  the  value  of  depended 

on  the  molecular  weight  of  the  original  polymer.  And  Clferri  and  Flory  [25] 

noted  that  in  the  case  of  natural  rubber,  a  specimen  croaslinked  by  gamma 

radiation  to  give  the  same  value  Cj  -  2.  0  kg/cm  as  one  crosslinked  with 

sulfur  with  a  2C^  =  2.  0,  exhibited  a  tower  value  of  2C^, namely,  l.  5  kg/cm^. 

Consequently  Krigbaum  and  Roe  [l8]urged  the  study  of  the  effect  of 
structural  differences  on  C^.  The  distribution  of  network  chain  lengths  could 
for  example  be  of  some  importance.  A  knowledge  of  the  factors  effecting  C^ 
would  help  in  correlating  experimental  results  such  as  the  energy  component 
of  the  elastic  force  obtained  by  different  workers  from  apparently  identical 
materials  [31,  32]  . 

Few  studies  have  been  carried  out  on  the  temperature  dependence  of  C^. 

When  C£  was  evaluated  from  stress -strain  data  at  very  short  times  apparent 
values  of  C2  decreased  with  increasing  temperature,  reflecting  a  faster 
approach  to  equilibrium  at  higher  temperatures  [28].  Under  conditions  where 
relaxation  effects  could  be  neglected,  the  temperature  coefficients  of  were 
found  to  be  positive  for  natural  rubber  [8]  ,  silicone  rubber  [28]  and  nearly 
zero  for  butyl  rubber  [28]  .  Roe  and  Krigbaum  [8,  30]  analyzed  stress -strain 
data  with  respect  to  entropy  and  internal  energy  contributions  to  elastic  force 
and  found  that  in  the  case  of  natural  rubber  it  is  the  internal  energy  effects 
which  are  largely  responsible  for  the  C£  term. 

Blatz  and  Ko[29l  working  with  polyurethane  foam  rubber  containing  50% 
air  by  volume,  found  that  2Cj  was  practically  negligible  and  that  all  of  the 
elastic  force  arose  from  the  2Cg  term.  This  finding  is  a  good  reminder  that 
the  existence  of  the  term  is  real  and  that  macroscopic  as  well  as  microscopic 
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structural  features  will  determine  the  relative  importance  of  the  C2  term. 

In  view  of  the  oversimplifications  inherent  in  the  statistical  theory,  one  can 
assume  that  all  the  structural  features  neglected  there  will  come  into  play  such 
as  volume  effects,  local  chain  packing  or  ordering,  non -Gaussian  segment 
distribution,  free  chain  ends,  crosslink  clustering,  non  equilibrium  and 
thermal  effects,  and  others.  In  view  of  these  complexities,  and  of  the  con¬ 
tradictions  so  far  encountered,  an  apriori  deduction  of  seems  rather  hope¬ 
less  at  this  stage,  and  that  efforts  should  be  concentrated  on  careful  measure¬ 
ments  of  values  of  as  a  function  of  rubber  structure. 
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F  Introduction  to  Dilatometry 

Our  knowledge  of  rubber  behavior  under  simple  elongation  would 
not  be  complete  without  an  insight  into  the  nature  of  the  changes  in  volume 
as  a  function  of  elongation  and  time.  From  a  theoretical  point  of  view,  the 
exact  function  of  allows  one  to  calculate  directly  the  quantites 

of  thermodynamic  interest  without  having  to  assume  certain  approximations 
as  has  been  the  case  in  previous  studies.  Following  the  change  in  volume 
at  a  fixed  strain  level  as  a  function  of  time,  should  permit  also  the  study  of 
ordering  or  crystallisation  kinetics,  since  at  higher  levels  of  order  the 
density  of  the  overall  sample  must  increase.  When  volume  changes  have 
been  used  to  follow  the  crystallisation  quantitatively  over  a  range  of  temp¬ 
erature  [331  ,  the  final  decrease  in  volume  on  crystallisation  was  found  to 
Ue  between  2.  0  and  2.  7  percent.  Experimentally,  measurement  of  A  V 
accompanying  elongation  was  first  attempted  by  Holt  and  McPherson  [  16*1 , 
who  found  no  volume  change  for  elongations  up  to  200%.  These  authors  were 
specifically  interested  in  effects  of  crystallisation  and  also  worked  with 
filled  samples. 

At  elongations  higher  than  200%  they  found  that  even  after  3  or  4  weeks, 
the  volume  of  the  stretched  rubber  decreased  at  an  approximately  uniform 
rate  with  the  logarithm  of  time.  These  decreases  were  the  greater,  the 
higher  the  elongation,  the  lower  the  temperature  or  the  longer  the  time 
the  rubber  was  kept  stretched.  These  observations  were  obviously  related 
to  the  kinetics  of  crystallisation.  In  later  work  [341,  Gent  followed  the 
crystallisation  of  the  sample  by  measuring  decreases  in  volume  and  re¬ 
laxation  of  the  stress.  He  found  that  the  degree  of  crystallisation  as 
measured  by  the  dilatometer,  was  proportional  to  the  reduction  in  stress. 

In  many  cases  the  equilibrium  value  of  stress  reached  was  sero,  in  some 
cases  the  test  piece  eventually  extended,  with  the  additional  extension 
amounting  to  up  to  5%  of  the  unstretched  length.  We  have  also  observed 
this  phenomena  of  relaxation  by  crystallisation,  but  only  in  the  lower  temp¬ 
eratures.  For  the  purpose  of  the  thermodynamic  studies,  our  lowest  temp¬ 
erature  was  30°C  and  the  highest  elongation  about  200%,  so  that  we  were  not 
in  the  crystallizing  region  of  elongation  and  temperature.  Lately  Wolstenholme 
measured  the  changes  in  volume  upon  elongation  of  a  series  of  cured  gum 
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Fig.  4,  ref  [35"'  :  Volume  change  measured  by  Wolstenholme 


as  a  function  of  elongation 
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elastomers  1*351  .  Using  a  water -fitted  dilatometer  submerged  in  a  thermo- 
stated  water  bath^naintained  at?  0.  015  °C,  Wolstenholme  claims  that 
changes  in  votume  of  10“®  cc  can  readily  be  determined.  He  notices  that 
the  votume  of  the  sample  changes  with  time  after  a  change  in  length,  which 
changes  were  characterized  by  a  rapid  drift  of  the  volume  in  the  first  half 
minute  after  stretching,  followed  by  a  very  slow  change  over  a  longer  time. 
Reproducible  results  were  obtainable  when  successive  elongation  changes 
were  made  every  three  minutes.  However,  water  as  a  confining  liquid, 
presented  the  additional  problem  of  swelling  of  the  sample  which  gave  rise 
to  a  slow  transient;  these  difficulties  were  minimized  by  pre-soaking  the 
samples  in  distilled  water  at  test  temperatures  for  a  minimum  of  24  hours. 
This  procedure  is  claimed  to  decrease  subsequent  water  absorption  by  the 
specimen  end  stabilizes  the  initial  sample  volume.  Subtracting  the  volume 
expansion  due  to  the  compressibility  of  rubber  from  the  total  volume  change, 
defines  an  induced  volume  change  which  Wolstenholme  attributes  to  the 
crystallization  of  the  sample.  He  tested  several  commerical  elastomer 
gums,  such  as  Hevea,  Neoprene,  Butyl,  SBR,  Paracril  B  and  Paracril  C. 
Elongations  applied  ranged  from  0  to  600%,  however  no  precise  data  was 
given  for  the  low  elongations.  Figure  4  shows  the  different  volume  changes 
of  a  Hevea  ear  pie  ■■  a  function  of  stressed  length.  The  compressibility 
correction  is  the  expansion  of  the  rubber  under  a  given  stress,  or  in  other 
words  defines  a  Poisson  ration,  p  smaller  than  the  incompressible  value 
of  0.  5.  In  filled  san  ss,  he  could  observe  an  apparent  Poission  ratio 
greater  than  0.  5  due  to  vacuole  formation  [  161.  Reduction  of  the  temperature 
from  35°C  to  0°C  lowers  the  elongation  at  which  crystallinity  becomes 
apparent  t35l.  Finally,  cycled  dilatometer  tests  lead  to  a  volume  hysteresis 
that  was  greatest  in  crystallizable  elastomers  and  practically  zero  in 
non-crystatlizable  elastomers. 

The  first  people  to  measure  the  volume  changes  in  detail  *t  low  and 
moderate  strains  were  Gee,  Stern  and  Treloar  [36]  .  From  thu  earlier 
data  of  Holt  and  McPherson  it  was  obvious  that  if  there  were  any  change  in 
volume  upon  elongating  the  sample  up  to  200%,  these  changes  would  have  to 
be  less  than  0.  2%  expressed  in  AV/V0  ,  and  Gee  [371  had  predicted  small 
increases  in  volume  at  low  strains.  For  any  isotropic  compressible  body  the 


tensile  force  mey  be  resolved  into  sheer  stresses  end  en  isotropic 
tension,  for  e  chenge  from  L  to  L  end  e  compressibility  K  ,  the  volume 
expension  AV  is  given  by  ^371*: 


L  (  —) 


dL 


P,  T 


(49) 


This  equetion  epplies  so  long  es  the  meterlel  remeins  isotropic,  end  is 
therefore  probebly  not  seriously  incorrect  for  rubber  et  elongetions  of  up 
to  100%.  By  meens  of  Equetion  (49).  it  is  thus  possible  to  predict  velues 
of  /V  from  the  stress -strein  curve,  or  else  one  cen  meesure  the  chenge  in 
volume  directly.  Insteed  of  e  diletometer,  Gee  end  co-workers  used  the 
method  of  hydrostetic  weighing  in  weter.  An  eccurecy  of  up  to  1  pert  in 
60,  000  in  volume  is  cleimed.  A  pure  gum  rubber  conteining  2%  sine  oxide 
geve  volume  chenges  epproximetely  double  of  those  expected  from  epplicetlon 
of  Equetion  (49).  This  devietion  wes  expleined  by  postuletlng  the  forming  of 
vecuoles  eround  the  filler  perticles,  which  merkedly  increese  the  volume  of 
the  semple.  Leter  test*?  were  run  on  e  semplo  crosslinked  by  meens  of  di- 
t -butyl  peroxide  end  showed  better  the  egreement  of  the  results  with 

the  theoreticel  predictions.  Figure  5  reproduces  Gee  end  co-workers  dete, 
plotted  egeinst  the  predicted  curve. 

Severet  yeers  leter,  Hewitt  end  Anthony  [  38  3,  redetermined  the  chenge 
in  volume  of  neturel  rubber  crosstinked  by  di-t-butyl  peroxide.  The  streins 
et  which  they  were  working  were  lower  then  those  used  by  Gee  end  his 
group  (40%  to  110%).  Agein,  to  obtein  high  sensitivity,  the  experimentel 
technique  wes  thet  of  meesuring  the  buoyency  in  weter.  The  tempereture 
control  wes  of  the  order  of  0.  005°C  es  meesured  with  e  Beckmenn  differentiet 
thermometer.  Figure  6  shows  the  percent  expension  for  streins  up  to 
^  s  l.  6.  Appllcetion  of  Equetion  (49)  results  in  the  theoreticel  curve  es 
seen  in  Fig.  6.  To  eveluete  (?)  f/3L)p,T  these  euthors  used  the  1st  term 
of  the  Mooney-Rivlin  equetion: 


o  27.5  “  32 ,5°C 
+  30  °C 

V  35  "  40#C 
o  45  “  50°C 
x  55“60*C 


^ig,  7,  ref  T9l  :  Internal  energy  contribution  from  Alien,  Bianchi  and  Price 
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and  integrated  (M/9L)  obtained  from  actual  stress-strain  measure- 

P>  i' 

ments.  Generally,  the  agreement  between  experimental  and  theoretical 
results  is  seen  to  be  good  . 

Allen,  Bianchi  and  Price,  r9]measured  the  change  in  volume  indirectly 
by  following  L  and  using  one  of  Maxwell *3  relationships: 

(3  f/?>  p )  :  (3V/3L)  (51) 

t,l  T,  p 

Experimentally  the  problem  consists  of  measuring  (d  f/dT)  at 
constant  volume,  that  is  to  say  necessitating  a  hydrostatic  ^  pressure 
of  some  120  atmospheres  over  a  range  of  5-10°  C.  An  accurate  estimation 
of  (df/dT)  was  obtained  by  the  increase  in  stress  at  approximately  1 

degree  intervals  over  5  degrees  in  all.  Indirect  measurements  of  (df/c:T) 

V,  L 

were  made  by  measuring  (a  f/Sp)^,  ^  ,  (dp/dTJy  L  and  (2tf/dT)p  L  . 


setting: 


fe  =  f  -  T(af/»T)LtV  (49) 

These  authors  obtain  f  If  as  a  function  of  elongation  (Fig.  7).  There  is 
by  no  means  good  agreement  with  the  results  of  other  workers  who  measured 
fe/f  indirectly.  In  any  case  these  results  confirmed  that,  at  moderate 
extension  ratios,  the  energetic  contribution  to  the  elastic  force  is  approximately 
20%  .  The  tests  were  not  at  sufficiently  low  stresses  to  confirm  Roe  and 
Krigbaum's  claim  that  contrary  to  the  behavior  of  a  Gaussian  network, 
there  is  an  energetic  contribution  to  the  elastic  force  which  varies  with 
extension  ratio  at  low  values  of  elongation. 

Finally  in  that  very  broad  approach  to  all  aspects  of  the  thermodynamics 
of  rubber-like  elasticity,  Allen  et  al.  [9l  calculate  the  dilations  of  their  samples 
as  a  function  of  elongation,  and  compare  the  results  with  their  data,  (Fig.  8). 

As  can  be  seen,  the  agreement  is  not  good,  but  never-the-less  the  authors 
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suggest  that  their  indirect  method  of  measuring  via  (<Jf/fcT)  would 
give  better  results  than  the  direct  measurement  of  dilationi  since  the  latter 
is  so  small. 

Concerning  the  use  of  water  as  a  liquid,  we  found  in  preliminary  studies 
that  the  weight  of  the  samples  Immersed  in  water  increased  markedly  as 
a  function  of  time.  Extrapolation  to  sero  time  showed  that  there  was  an 
immediate  gain  in  weight  upon  immersion.  The  use  of  swelling  or  inter¬ 
acting  fluids  is  of  course  highly  questionable  because  of  the  introduction 
of  new  complexities;  thus  the  choice  of  a  non  interacting  fluid  will  be  discussed 
in  the  experimental  section. 

Rough  estimates  even  suffice  to  show  that  the  dilatometer  cavity  acts 
as  an  extremely  sensitive  thermometer.  With  a  capillary  of  0.  04  cm 
diameter,  it  is  possible  to  detect  changes  of  0.  0001°C  using  an  accurate 
cathetometer.  In  view  of  the  work  by  Wolstenholme,  using  water  as  confining 
liquid,  we  thought  first  to  set  the  bath  temperature  at  3.  97°C  (maximum 
density  of  H^O)  so  that,  if  water  were  used  as  the  dilatometer  fluid  since  its 
coefficient  of  thermal  expansion  is  zero  at  this  temperature  there  should  be 
no  thermometer  effect  on  the  dilatometer.  Since,  subsequent  work  showed 
that  water  should  not  be  used  as  the  confining  fluid,  this  line  had  to  be 
abandoned. 
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ft,  "Fast"  Stretching 

The  temperature  rise  during  adiabatic  etretching  of  rubber  has  been 
known  for  a  long  time.  The  first  observations  were  reported  by  Gough 
in  1805  fl]  ,  but  it  was  Kelvin  who  first  considered  the  thermodynamic 
implications  and  suceeded  in  showing  that  the  positive  stress -temperature 
coefficient  was  a  necessary  thermodynamic  consequence  of  the  evolution 
of  heat  on  extension.  A  simple  application  of  the  second  taw  of  thermodynamics 
to  rubber  deformation  teads,  as  stated  earlier,  to 

■“W  =  T<1S  (52) 

Thus,  since  heat  is  given  off,  and  the  sign  of  d&is  negative,  the  entropy 
must  decrease  presumably  due  to  the  aligning  of  the  chains  and  thereby  there  is 
toss  of  randomness. 

There  are  three  different  ways  of  measuring  the  heat  given  off  by  a 
rubber  during  stretching. 

(a)  direct  calorimetric  measurement  of  the  heat  exchanged  when 
rubber  is  stretched  isothermally;  (b)  determination  of  the  temperature 
rise  during  "adiabatic"  stretching,  and  (c)  measurement  of  the  change 
in  heat  content  between  stretched  and  unstretched  rubber. 

All  work  to  date  has  been  done  by  method  (b).  In  the  method  used  by 
Dart,  Anthony  and  Guth  [39l  ,  two  identical  rubber  samples  were  held  to¬ 
gether  by  twisting  the  samples  and  a  thermocouple  bead  placed  between  them. 
The  authors  found  that  usually  the  temperature  change  on  extension  was 
different  from  that  on  retraction.  Generally,  samples  which  were  capable 
of  crystallizing  exhibited  a  marked  temperature  rise  at  the  onset  of  crystal¬ 
linity  and  there  was  even  less  reversibility  during  the  stretching- 
relaxing  process,  and  the  cooling  on  retraction  was  generally  greater  than 
the  heating  on  extension.  The  last  phenomenon  was  attributed  to  a  time  lag 
in  the  crystallization  process  since  as  the  time  during  which  the  rubbers 
were  kept  stretched  before  retraction  was  increased  the  differences  between 
the  two  heats  become  larger.  This  was  taken  to  show  that  crystallization 
took  place  before  retraction.  In  this  context  it  is  of  interest  to  note  that  the 
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Fig.  9,  Ref  fl2l  :  Temperature  rise  in  adiabatic  extension  from  Treloar 


Ill  -  36 


temperature  response  was  imznadlata  for  slongations  of  last  than  300  par 
cant,  but  at  largar  onas  thara  was  a  lag  of  tha  tamparatura  [39l  which 
amountad  to  as  much  as  5  saconds  for  high  alongations  and  could  not  ba  dua 
to  slowness  in  tha  galvonomatar  of  thermocoupla  rasponsa.  Tha  authors 
stats:  "This  impliss  soma  procass  which  takas  a  small  amount  of  time  must 
ba  going  on  in  tha  rubbar  during  elongation".  One  might  suppose  that  there 
is  a  diffusion  process  by  which  the  chains  accomodate  themselves  to  the 
strained  positions. 

Figure  9  shows  the  temperature  rise  in  so-called  "adiabatic"  extension 
compiled  from  the  earlier  measurements  undertaken  by  Joule,  James  and 
Cuth  T  12 ].  The  data  shows,  at  small  elongations,  an  initial  cooling  followed 
by  rapidly  rising  warming  affect  as  the  extension  was  increased.  This  course 
of  the  thermal  changes  must  be  compared  to  the  entropy  changes  derived  from 
stress -temperature  measurements.  Thus,  the  initial  cooling  corresponds 
to  the  initial  positive  entropy  of  extension  associated  with  the  expansivity  of 
the  rubbers  and  breaking  of  chain  associations  while  the  subsequent  heating 
corresponds  to  the  large  negative  entropy  term  from  the  drop  in  number  of 
available  conformations  of  the  network  molecules.  At  still  higher  elongations 
possible  crystallisation  phenomena  can  introduce  further  and  larger  heat 
effects.  It  is  of  interest  that  rubbers  which  can  not  crystalline  do  not  heat 
as  extensively  as  those  that  are  capable  of  crystallisation.  For  example,  one 
can  obtain  differences  in  temperatures  of  14°C  in  natural  rubber  at  600% 
elonation,  but  only  2°C  in  a  butadiene-acrylic  rubber  at  the  same  extension. 

In  this  connection,  fast  stretching  thermal  measurements,  which  are  less 
accurate  than  stress -strain  measurements,  point  out  in  a  very  direct  manner 
to  a  necessary  consequence  of  the  kinetic  thecrv  of  elasticity:  namely  that 
"the  process  of  deformation  of  rubber  is  capable  of  a  reversible  transformation 
of  work  into  heat"  [121 

For  "fast-stretching",  the  thermodynamics  wilt  be  derived  as  under 
"adiabatic".  In  the  Discussion  it  witl  be  shown  why  these  two  conditions 
are  hardly  strictly  equivalent.  Under  truly  adiabatic  conditions,  the  net 
entropy  change  molecular  ordering  plus  self-heating  is  zero,  and  the  change 
in  temperature  dT  is  given  by: 

dT  3  dH/C  .  (53) 

p, 
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where  dH  ia  tha  amount  of  haat  davalopad  in  tha  procaaa,  and  C  .  tha 

D  |  Ij 

•pacific  h«at  at  constant  lsngth, 

Tha  total  changa  in  intarnal  anargy  (adiabatic  and  at  conatant  praaaura)  ia 
given  by: 


/E 


tot. 


/  f  dL  -A  /  TdS  -  Jpd 


(54) 


which  can  ba  written  in  ita  differential  form 


*Etot,  :  t  +  T(JS_j  .  -  p(^V)  (55) 
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But  from  tha  Maxwell  relationahipa,  wa  have  aaan  that 
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If  tha  rubber  ia  deformed  adiabatieally,  than  ainca  dS  -  0  by  definition, 


/IE, 


fdL 


pdV 


(57) 


And  now  if  the  ayatem  ia  allowed  to  tranamit  the  heat  developed  during  tha 
adiabatic  atretching  proceaa,  then 


(58) 


But  the  pathway  a  which  tead  to  the  combined  energy  chaige  -f-  A E 

1  2 

leavea  the  aample  in  tha  aama  final  condition  aa  in  the  total  proceaa 
aaaociated  with  Hence  AE  4~  AE^  and 
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and  we  obtain  finally, 


(60) 


If  the  further  assumption  is  made,  that  Cy(  L 
range  in  question,  then 

f(  Sr' 

.  tJ  StPil 

*T  = _ 

c 

V ,  L 


is  independent  of  T  over  the 


(61) 


However,  as  in  most  cases  when  dealing  with  a  liquid  or  a  solid,  the 
difference  between  the  internal  energy  and  the  enthalpy,  or  the  work  term 
pdV,  is  considered  negligible.  Thus  we  can  write  equally  well: 


at 


{ (*i) 

TJ  3T  .  dL 
J  p,  L, 


(62) 


which  has  been  given  by  Leaderman  [40  ■  in  the  differential  form: 
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SECTION  I 


FORCE  VERSUS  TEMPERATURE 


MEASUREMENTS 


Fig.  11  Overall  Experimental  Layout  for  Force  versus  Temperature  and  "fast-stretching1 
experiments  for  Rubber. 
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X.  EXPERIMENTAL 

a.  igiiBua&RiiflgR  aad.  mtahinUil  dmelstf 

In  m  study  of  the  responses  of  rubbers  to  mechanical  deformation 
analysis  of  the  results  will  require  the  use  of  thermodynamic,  elastic  and 
rheological  equations  of  state.  Consequently  it  will  be  necessary  to  re¬ 
present  the  variables  in  3-dimensional  space  where  the  thermodynamic 
functions  as  ordinate  will  be  placed  along  the  s-axis,  whilst  the  x-and 
y-axis  will  show  the  termperature  and  the  extension  ratio.  Thus  the 
thermodynamic  state  functions  will  be  treated  as  dependent  variables. 

During  the  force  versus  temperature  measurements,  the  force 
is  directly  measured  by  a  calibrated  load  cell  or  more  accurately  described 
as  a  transducer.  The  instrument  used  in  all  the  tests  was  supplied  by 
STATHAM  Instruments,  Los  Angeles  64,  model  number  Gl-80-350.  The 
range  of  load  is  +-  80  os.  The  excitation  potential  is  15  volts,  and  the 
calibration  factor  is  given  as  53.  08  microvolts  (open-circuit)  per  volt 
per  os.  This  unit  is  not  temperature  compensated,  and  therefore  shows  a 
thermal  zero  shift  and  a  thermal  sensitivity  shift.  Since,  the  transducer 
should  be  kept  at  constant  temperature  during  the  experiments,  thermal 
insulation  of  the  transducer  was  necessary.  An  Invar  drop  rod  was  used, 
to  hook  up  the  sample  to  the  load  cell  (Fig.  10). 

The  strain  sensitive  wire  elements  of  the  transducer  are  arranged 
in  the  form  of  a  Wheatstone  bridge.  Direct  current  was  used  to  excite  the 
transducer.  The  value  of  the  potential  applied  must  be  a  constant  in  time. 

A  HARRISON  Laboratories,  Berkeley  Heights,  New  Jersey,  Model  801 
C  Power  Supply  was  used  in  conjunction  with  a  SOLA  Electric  Company, 
Chicago,  constant  voltage  transformer,  catalogue  number  30804.  This 
input  potential  circuit  proved  to  be  remarkably  stable  (of  the  order  of 
1  mv  out  of  15  v  over  a  period  of  3  months). 

The  force  versus  temperature  measurements  being  carried  out 
at  constant  length,  both  ends  of  the  sample  holders  are  clamped  down,  and 
only  the  force  varies  as  a  function  of  temperature.  The  sample  is  con¬ 
tained  in  a  large  converted  refrigerator  (Fig.  11  )  which  can  aither  be 
maintained  at  constant  air -temperature  at  any  point  between  -  25°C  and 
+  130°C,  or  can  serve  as  an  environmental  chamber  as  the  temperature  is 
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varied.  The  time-temperature  profile  is  not  linear,  but  as  a  complete 
temperature  cycle  is  of  the  order  of  12  hours,  thereby  thermal  equilibrium 
throughout  the  sampel  is  ensured  at  all  times.  An  added  feature  is  that 
the  atmosphere  of  the  chamber  can  also  be  controlled;  in  most  of  our  tests, 
nitrogen  gas  surrounded  the  sample  continuously. 

A  compressor  and  a  series  of  strip  heaters  allow  one  to  reach 
any  temperature  in  a  given  range.  Most  of  the  force  versus  temperature 
experiments  were  carried  out  overnight  in  a  time  cycle  of  temperatures 
from  2  5°  C  to  60°  C.  The  heaters  are  wired  to  large  15  ampere  power- 
stats  ,  and  it  is  possible  to  read  the  amperes  of  current  which  are 
supplied  to  the  various  heaters  directly.  It  is  also  possible  to  calibrate 
the  temperature  of  the  environmental  chamber  at  equilibrium  as  a 
function  of  the  power  supplied  to  the  heaters.  As  expected  (FigJ2  ),  these 
curves  are  roughly  linear,  since  the  heat  evolved  is  proportional  to  the 
power  which  is  in  turn  proportional  to  the  square  of  the  current.  These 
calibrations  are  extremely  useful  in  obtaining  the  desired  temperature, 
or  range  of  temperatures.  To  enable  nightime  testing,  a  spring  loaded 
relay  w.  s  installed  in  the  heating  circuit;  as  the  temperature  rises  above 
a  certain  predetermined  value,  the  relay  opens  the  heating  circuit,  and 
the  temperature  slowly  decreases  by  losses  to  the  ambient  as  a  function 
of  time.  Good  thermal  insulation  ensured  the  obtention  of  a  long  tempera¬ 
ture-time  cycle  and  thermal  equilibrium  within  the  samples. 

A  copper -constant an  thermocuple  from  the  OMEGA  Engineering, 

Box  47,  Springdale,  Connecticut  was  used  to  measure  the  temperatures. 

The  reference  temperature  of  0°  C  was  established  by  the  physical  equi¬ 
librium  of  ice  and  water,  sustained  by  a  thermoelectric  cooler  automatic¬ 
ally  controlled  by  a  bellows -micro3witch  sensing  mechanism  which  res¬ 
ponds  to  the  relative  volumes  of  ice  and  water  in  the  hermetic  cell.  This 
Ice-Point  Thermocouple  Reference  System. is  manufactured  by  the 
Joseph  KAYE  and  Co.,  Cambridge  38,  Massachusetts;  model  1150,  and 
must  be  specified  for  correct  use  of  the  1160  interchangeable  Thermocouple 
Probe  Assembly,  also  manufactured  by  the  Joseph  KAYE  Company. 

An  x-y  Potentiometric  Recorder,  of  the  HEWLETT-PACKARD, 
MOSELEY  division,  Pasadena,  California,  model  7001  AM  was  used  to 
determine  the  force -temperature  characteristics  of  the  rubber.  On 
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ig.  13  Sample  holder  and  detail  of  self -tightening  jaws 
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potentiometric  mode,  large  input  resistances  are  allowable  in  the  most 
sensitive  ranges.  Temperature  differences  of  the  order  of  0,001°  C, 
and  differences  in  stresses  of  the  order  of  1%  were  easily  analyzed.  The 
pen  system  used  was  a  Cartridge  loader,  which  after  some  adjustments, 
proved  to  be  satisfactory. 

A  short  description  of  the  clamping  system  follows.  Previous  very 
extensive  efforts  had  failed  to  devise  an  adhesive  bond  which  would  hold 
at  elevated  temperatures  and  elongations  t47^  Figurei3  shows  the  type  of 
self-tightening  end  clamps  that  were  finally  constructed.  The  rubber 
sample  is  imbedded  between  two  rubber  strips  tightly  held  by  the  self- 
locking  triangular  aluminum  wedges.  As  pressure  is  applied  to  the  wide 
surfaces  of  these  wedges  via  the  back-pressure  screw,  the  wedges  tend 
to  slide  forward  and  lock  the  sample  in  its  grip.  While  this  method  of 
l»ft|di«*g  the  sample  is  highly  satisfactory  at  moderate  and  high  levels  of 
strain,  the  inaccuracies  of  this  system  become  quite  large  at  very  low 
values  of  strain  (X.  <  1.  1,  where  X  is  the  extension  ratio).  This  fact  will 
be  illustrated  further  on,  in  the  discussion  of  the  weighting  factor  to  be 
applied  to  the  Mooney -Rivlin  type  plot  (See  Appendix  1  ). 

For  the  measurement  of  the  elongation  of  the  sample,  an  overall 
length  was  measured  between  the  pads  of  rubber.  The  rest  length  was 
estimated  in  the  same  manner.  Then  after  the  test,  the  rubber  was  taken 
out  and  a  mark  appeared  on  the  rubber  at  the  place  of  the  pads.  This 
total  length  between  pad  marks  was  taken  for  the  zero  length,  and  all 
other  lengths  were  corrected  by  the  difference  in  the  two  initial  lengths 
(these  were  not  appreciably  different).  A  trial  run  was  carried  out  by 
marking  6  lines  on  the  rubber  with  white  ink,  and  measuring  the  relative 
distances  between  the  markings  with  a  cathetometer .  Plotting  the  various 
positions  as  a  function  of  the  initial  positions  (with  the  rubber  at  rest), 
and  applying  a  least  square  linear  fit  to  the  points  so  obtained,  the  slopes 
of  these  lines  should  correspond  to  the  different  elongations,  measured 
by  taking  overall  lengths.  The  "best"  elongations  obtained  by  least 
squares  analysis  correspond  fairly  closely  to  the  elongations  obtained 
from  an  analysis  of  the  overall  lengths,  however  the  standard  deviations 
are  quite  high,  being  of  the  order  of  5  to  10%.  A  later  section  will  be 
devoted  to  the  problems  of  sample  shape  and  sample  holders. 
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ELONGATION  *  328% 
TEMPERATURE  *  57.5° 
No  ATMOSPHERE 
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B.  Experimental  Technique 
1,  Stress  relaxation 

Adopting  the  technique  of  previous  workers  [420  *  the  samples  were 
held  for  a  long  period  of  time  at  the  highest  temperature  and  at  the  highest 
elongation  to  insure  minimal  relaxation  of  the  stress  as  a  function  of  time 
during  the  actual  time -temperature  cycles  of  the  tests.  A  blanket  of  Ng 
gas  is  kept  surrounding  the  sample  at  all  times  to  prevent  oxidative  de¬ 
gradation  at  the  higher  temperatures. 

The  total  force  on  the  sample  is  then  read  from  the  xy  recorder  and 
plotted  as  a  function  of  the  logarithm  of  time.  There  is  no  claim  made 
as  to  the  mechanism  of  relaxation,  whether  it  be  one  of  the  forms  of  creep, 
crystallization  or  agings  However,  if  the  plot  of  F  =  f  (log^g  4  t)  where  F 
is  the  tctal  force  is  linear,  then  as  can  be  seen  from  Fig.  14  and  FigJL5, 
it  is  easy  to  calculate  the  time  after  which  testing  can  commence,  allowing 
a  two  percent  overall  relaxation  of  force  at  the  highest  temperature  and 
elongation  over  the  entire  period  of  the  tests. 

From  Fig  J4  ,  the  equation  of  the  relaxing  specimen  is  given  by: 


F  =  A  log1Q  At  +  B 


(64) 


where  F  is  the  total  force,  and  A  and  B  are  the  linear  coefficients.  The 
coefficients  A  and  B  can  be  read  from  the  plot.  Let  us  now  solve  for  tj, 
where  At  =  t^  -  t^  =  amount  of  time  for  total  test,  and  where  cF  =  (F2-Fj)/F^ 
s  0.02  :  cF  is  the  relative  error  on  F  .  From  (64) 

A(l°gl0  (\  +At)  -  lo810  V 


cF  = 


(65) 


A  log10  tj  +  B 


Rearranging  (65),  log1Q  tj  =  log^l  +At* 


BeF 


1  -  c  F 


A(l-eF) 


(66) 


This  equation  can  be  solved  graphically,  as  shown  in  Fig.15  .  Using 
this  technique,  it  is  very  easy  to  know  when  to  Btart  the  actual  testing. 
Once  the  nominal  relaxation  for  the  period  of  testing,  is  less  than  two 
percent  at  the  highest  elongation  and  temperature,  then  the  relaxation  is 
considered  sufficiently  accomplished,  and  the  temperature  is  decreased 
by  shutting  off  the  heaters,  and  the  sample  is  allowed  to  equilibrate  at 
room  temperature  and  in  the  environment  of  N2.  Using  the  previously 
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described  technique  of  heating  and  then  allowing  the  sample  to  cool  over 
a  period  of  approximately  12  hours,  a  plot  of  F(force)  expressed  in  milli¬ 
volts  versus  T(temperature)  also  in  millivolts  is  obtained  on  the  xy 
recorder,  as  a  series  of  F  versus  T  lines  at  constant  length  of  the  sample. 
Using  thermocouple  calibration  constants  given  in  the  literature[43]  ,  the 
temperature  can  be  expressed  in  degrees  Centigrade  or  Absolute.  The 
calibration  of  the  load  cell  is  done  by  hanging  known  weights  on  the  sample 
holder,  and  plotting  the  reduced  output  (output  divided  by  potential  applied) 
as  a  function  of  the  known  weights.  Plotting  the  apparent  weights  at  first 
allows  one  to  extrapolate  to  the  dead  weight  of  the  holder  and  Invar  drop 
rod.  The  second  line  of  Fig.  16  shows  the  true  calibration  curve.  A 
digital  voltmeter  was  used  to  check  from  time  to  time  the  input  potential 
to  the  load  cell.  It  was  found  to  be  constant  to  lmv  (out  of  15  v  applied) 
over  a  3  month  period  of  time. 

Knowing  all  of  the  calibration  constants  (thermocouple  and  load  cell), 
and  the  zero  values  and  ranges  of  the  x  and  y  axes,  a  simple  computer 
program  (P  31)  calculates,  for  each  length,  the  forces  (dynes/cm2)  and 
temperatures  (degrees  absolute,  degrees  centigrade).  A  listing  of  the 
Fortran  IV  source  deck  as  well  as  some  typical  entries  and  outputs  (both 
written  and  punched)  are  discussed  in  Appendix  I.. 

At  each  length,  two  time -temperature  cycles  were  carried  out  to 
check  upon  the  reversibility  and  reproducibility  of  the  data.  Figjl7  shows 
F  versus  T  for  four  different  elongations.  It  can  be  seen  that  the  F-T 
diagrams  are  reproducible,  and  except  at  the  very  highest  elongations 
(350%),  all  runs  were  completely  eversible.  In  order  to  save  time,  sub¬ 
sequent  testing  would  only  be  carried  out  on  a  single  cycle  per  length. 

From  this  point  on,  all  of  the  raw  data  has  been  collected  and  now 
only  needs  to  be  processed.  A  series  of  consecutive  computer  programs 
has  been  written  to  calculate  the  various  quantities  of  interest.  Each 
punched  output  is,  with  some  slight  modifications,  the  input  for  the  next 
program.  Thus  a  complete  analysis  of  the  various  parameters  at  any 
stage  of  the  study  is  easy. 

From  the  force-temperature  data  at  constant  length,  it  is  possible 
to  calculate  the  stress -strain  curves  at  constant  temperature.  But  these 
curves  must  be  corrected  for  the  thermal  expansion  of  the  rubber,  thus 
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obtaining  corrected  stress-strain  data.  It  is  possible  to  then  recalculate 
corrected  force -temperature  data  at  constant  length.  One  can  also  fit 
the  data  to  a  Mooney-Rivlin  type  plot.  It  is  all  of  these  calculations  that 
are  done  in  the  successive  computer  programs. 

2.  Sample  preparation 

The  samples  under  discussion  all  came  from  the  same  supply  (1/ 16" ) 
of  DUNLOP  NATURAL  RUBBER  with  one  percent  dicumyl  peroxide  as 
the  crosslinking  agent  and  were  obtained  by  courtesy  of  Dr.  DINGLE.  The 
samples  were  stored  cold  in  the  dark  and,  as  much  as  }>>ssible,  all 
experiments  were  carried  out  on  specimens  from  the  same  slab. 

The  actual  geometrical  configuration  of  the  specimens  used  were 
rectangular  pieces  80  mm  x  2  4  mm  cut  with  a  die  manufactured  by  the 
CLEAR  CUT  STEEL  RULE  DIE  Corp.,  Plainview,  Long  Island.  Great 
care  had  to  be  taken  in  the  cutting  of  the  samples,  in  order  to  avoid  nicks 
or  tears  which  serve  as  local  stress  concentration  points  and  initiate 
rupture  of  the  sample.  A  small  hydraulic  press  was  used  in  conjunction 
with  a  slab  of  a  soft  plastic  placed  between  the  press  and  the  sample  to 
be  cut.  The  cut  sample  was  then  carefully  measured  at  different  points 
with  a  micrometer  from  STARRETT  Number  1010M  as  to  uniformity  of 
thickness  and  width.  The  standard  deviations  on  thickness  and  width 
were  of  the  order  of  three  percent.  The  length  of  the  sample  before 
placing  it  in  its  grips  was  measured  with  a  ruler  calibrated  in  millimeters. 
A  reading  could  be  easily  estimated  to  within  0.  5  mm. 
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FIG.  18  RELAXATION  CURVE  FOR  OAI-22 
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TABLE  I 


Number  of 

t 

lo*10  * 

net  F 

Point 

(hra) 

(cms) 

0 

0 

m 

• 

1 

0.33 

-0.478 

23.80 

2 

0.66 

-0.  181 

23.  55 

3 

1 

0 

23.45 

4 

1.50 

0.  176 

23.40 

5 

4.33 

0.636 

23.75 

6 

21.  12 

11.324 

22.85 

7 

43.62 

1.639 

22.65 

Stress  Relaxation  of  DA  1-22  as  a  function  of 

of  time 

TABLE  II 

Number  of 

t 

lo«10  * 

net  F 

Point 

(hrs) 

(cms) 

0 

0 

- 

1 

0.25 

-0.699 

23.40 

2 

0.50 

-0.301 

23.30 

3 

0.75 

-0. 125 

23.20 

4 

1.00 

0 

23.  15 

5 

15.50 

1.  190 

23.70 

6 

19.33 

11.286 

22.65 

7 

49.25 

1.692 

22.30 

Stress  Relaxation  of  DA  1-23  aa  a  function 


of  time 
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FIG.  20  RELAXATION  CURVE  FOR  DA- 15 
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FIG.  21  RELAXATION  CURVE  FOR  DB- 
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FIG.  Z2  RELAXATION  CURVE  FOR 
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II  RESULTS  - 

Two  samples  of  natural  rubber  were  analyzed  exhaustively  and  these 
results  shall  now  be  presented.  The  samples  examined  are  numbered 
DA  1-22  and  DA  1-23.  All  measurements  taken  after  relaxation  and  con¬ 
secutive  testing  cycles.  The  dimensions  are  given  below: 


Rest  length 

Thickness 

Width 

Area 

(mm) 

*mm) 

(mm) 

(cm2) 

DA  1-22 

79 

1.89 

4,11 

7.603  .10" 

DA  1-23 

77 

1.82 

4.06 

7.  389  .10 

B.  Relaxation  of  Samples 

The  samples  were  each  held  at  about  175%  elongation  over  a  period 
of  several  days,  at  the  highest  temperature  of  50°C.  The  force  or  total 
stress  on  the  sample  was  noted  as  a  function  of  log^g  t,  where  t  is  the 
time  in  hours.  From  Figures  18,  19  and  Tables  I,  XI,  it  is  evident 
that  the  relaxation  of  stress  is  linear  with  log^g  t  .  It  was  found  on  pre¬ 
vious  samples  that  either  at  higher  temperatures  (T  >  50°C),  or  at  higher 
elongations  (\  >3.00),  the  relaxation  curves  were  initially  linear  but 
then  the  stresses  sharpl’  diminished  with  increasing  times  (Fig20  ). 
Sometimes  the  samples  even  broke  during  the  relaxation  tests.  Higher 
temperatures  would  naturally  increase  relaxation,  whether  it  be  creep 
due  to  viscoelastic  or  to  chemical  relaxation  [  441  .  Higher  elongations 
would  of  course  increase  relaxation  induced  by  crystallization,  as  well 
as  accelerating  any  crack  propagation.  We  were  notably  unable  to 
"relax"  the  synthetic  natural  rubbers  (NATSYN)  or  the  Shell  elastomer 
(Figs.2l,22)due  to  failure  after  approximately  100  hours.  However, 
probable  linear  relaxation  could  be  achieved  by  working  at  substantially 
lower  elongations  and  lower  temperatures  with  both  of  these  elastomers. 

C.  Treatment  of  data 

In  this  section  a  brief  description  of  the  various  calculations  carried 
out  upon  the  raw  data  by  means  of  a  computer  program  will  be  given. 

From  there  a  series  of  consecutive  programs  carry  out  the  prescribed 
algebraic  operations.  From  the  F  versus  T  data,  as  mentioned  in  the 
introduction,  it  is  possible  to  calculate  various  changes  in  state 
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function.  However,  it  is  also  possible  to  present  the  data  in  the  form 
prescribed  by  MOONEY  and  RIVUN  [45}. 

It  must  be  remarked  at  the  onset  that  all  of  the  tests  are  carried 
out  at  constant  length.  The  thermodynamic  relationships  are  a  function 
of  ratio  of  lengths  at  a  given  temperature.  At  elevated  temperatures, 
the  normalized  length  is  the  total  length  divided  by  the  rest  length  at 
that  same  temperature.  If  the  slope  of  the  force  versus  temperature 
data  were  to  be  plotted  as  a  function  of  length,  one  would  find  at  very 
low  elongations  negative  values  of  the  slope  ( Fig,  23).  This  is  the  thermo  - 
elastic  inversion  phenomena,  first  pointed  out  by  Anthony,  Caston  and 
Guth  [46]  •  It  is  easy  to  calculate  a  corrected  elongation  on  the  basis  of 
the  unstressed  length  at  the  temperature  in  question  rather  than  at  room 
temperature.  A  set  of  stress-strain  values  is  obtained  at  constant 
length.  Each  value  of  strain  is  then  corrected  for  thermal  volume 
expansion  by  means  of  the  equation; 

Xc  =  X  e“  *  A  T  (67) 

where  a  is  the  mean  linear  temperature  coefficient  of  the  unstrained 
rubber  and  A  T  is  the  temperature  rise  in  °C  above  room  temperature. 
Equation  (j  )  Is  the  most  general  form  expressing  the  temperature  de¬ 
pendence  of  a  given  parameter.  If  <T  were  to  be  replaced  by  a  =  a(T), 
then  in  order  to  calculate  the  temperature  dependence  of  the  parameter, 
one  would  have  to  integrate  Equation(  $7)overall  values  of  T. 

Recently  Allen,  Bianchi  and  Price [9  ]  showed  that  a  is  independent 
of  elongation  and  temperature  in  the  intervals  \  =  1.0  -2.2  and  T=  3<f  C 
-70°C.  We  shall  use  their  value  namely  a  =  2.196  x  10~^  deg”**.  This 
is  the  only  correction  that  these  or  other  worker"5  in  the  field  mention. 
However  as  pointed  out  by  Blatzf47]  there  is  also  a  correction  fer  the 
stresses.  These  are  forces  per  unit  area,  and  the  area  is  equally 
temperature  dependent.  Hence  this  work  comprises  also  correction 
terms  on  the  stresses:  namely, 

-2oAT 

%  =  (68) 

However,  these  corrections  are  very  small  and  most  authors  simply 
neglect  them  [  42] 
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In  essence,  the  computer  programs  were  written  in  order  to 
assemble  and  index  the  various  data,  to  effect  cuts  in  the  implied 
3-dimensional  surface  of  F-T-X  perpendicular  to  the  axis,  to  correct 
for  thermal  expansion  following  Equations  16?)  and  (68)  the  stress-strain 
data,  and  finally  to  recompute  on  the  one  hand  the  corrected  force - 
temperature  data  from  which  the  changes  in  state  functions  may  be 
derived,  and  on  the  other  hand  to  calculate  the  data  in  form  of  the 
Mooney-Rivlin  typeplot.  These  various  programs  are  described  in 
Appendix  I,  where  a  general  summary  of  the  programs  used  is  given, 
as  well  as  a  detailed  analysis  of  each  one. 

It  is  interesting  to  note  that  the  least  square  analysis  is  &  weighted 
least  square,  where  the  weight  ascribed  to  a  point  is  the  reciprocal  of 
the  variance  of  that  point.  But  onr.;>  th •*  original  data  are  considered  to 
be  of  equal  weight,  all  of  the  subsequent  weights  are  defined.  This 
allows  for  much  more  precise  interpretation  of  the  results.  In  the  least 
squares  fitting  technique,  we  assume  no  error  in  the  independent  variables. 
This  is  not  an  essential  assumption,  but  simplifies  the  computations 
greatly. 

In  Appendix  I  ,  the  explanation,  listing  and  input-output  formats 
will  be  described.  Thus  the  actual  results  will  be  presented  in  tabular 
form  just  as  they  appear  on  the  computer  output  sheets 
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HI,  DISCUSSION 

Foret  viraui  Temperature 

Our  study  of  the  force-ver sus-temperature  data  for  tha  two  specimens 
analysed  brings  out  a  numbar  of  intaraatlng  points.  Tha  first  phanomenon  of 
intarast  ara  some  aspects  of  tha  relaxation  curves  (Sea  figs  18,  19).  As  already 
mentioned  in  the  Introduction,  [26,  27,  28)  where  as  the  classical  viscoelastic 
relaxation  la  given  by  a  linear  plot  of  log  force  as  a  function  of  time,  we  find 
here  a  linear  plot  of  the  decay  of  force  as  a  function  of  logarithm  time.  The 
same  was  also  found  to  be  tha  case  in  a  study  of  tha  thermodynamic  interpretations 
of  tue  elastic  properties  of  rubbers  obtained  from  Ethylene-Propylene  copolymers 
[493  .  Why  the  relaxation  should  obey  this  law  is  not  clear;  a  further  study  would 
be  of  great  interest.  It  is  well  established  however,  that  in  order  to  ensure 
reversibility  of  the  stress  temperature  runs,  the  force  at  a  given  elongation  and 
temperature  must  remain  constant  to  within  the  tolerances  of  our  experimental 
technique.  Few  authors  specify  the  cut  off  point  where  the  relaxation  becomes 
too  large  so  that  we  endeavoured  to  set  a  point  arbitrarily.  Thus  we  took  a  total 
decrease  in  force  over  the  entire  period  of  testing,  of  less  than  1  %( at  the  highest 
level  of  temperature  and  elongation)  as  applicable.  A  preliminary  run  of  force 
temperature,  carried  out  after  such  a  relaxation,  showed  that  sucessive  runs 
were  undistinguishable.  However  since  the  second  term  of  the  Mooney-Rivlin 
equation  (2CZ)  is  felt  by  some  authors  to  arise  from  th*  lack  of  attaining  equil¬ 
ibrium,  by  this  criterion  however  small  the  relaxation,  this  procedure  could  not 
be  used  to  ascertain  the  existence  of  C ^  at  true  equilibrium.  Before  a  decision 
is  possible,  one  would  have  to  wait  for  a  molecular  understanding  of  the  relaxation 
phenomenon  to  give  us  better  Insight  into  the  true  meaning  of  the  2C^  term,  if 
such  a  term  exists  in  a  true  equilibrium  state. 

In  Appendix  I  it  has  been  shown  how  very  sensitive  the  function  rr/X  -  X  is 

when  applied  to  the  Mooney-Rivlin  plot.  Van  der  Hoff  [22"'  indeed  corroborates 

our  findings  doing  the  error  analysis  of  our  and  his  data,  by  showing  the  effect 

of  a  very  small  change  in  rest  length  upon  any  fitting  to  the  Mooney-Rivlin  plot 

(Fig.  24).  It  can  be  clearly  seen  that  if  L.  is  not  known  with  great  accuracy, 

o 

the  approach  of  the  Mooney-Rivlin  plot  to  X*4  1,  is  meaningless.  Van  der  Hoff 
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asserts  that  extremely  careful  measurements  of  streaa  versua  strain,  and  a 
least  squares  linear  extrapolation  to  find  the  rest  length,  allows  one  to  measure 
a  value  of  2C,  -  0  whan  X  is  close  to  unity.  We  had  found  a  substantial , 

increasing, departure  from  linearity  close  to  the  rest  elongation  for  both  of  our 
samples.  This  departure  was  taken  to  mean  thrt  the  rest  lengths  were  not 
accurately  known.  Consequently  in  the  sensitive  range  of  the  Mooney-Rivlin 
plot,  the  data  were  effectively  discarded  by  the  use  of  a  weighting  function  as 
described  in  Appendix  I. 

Similarity  a  careful  analysis  of  the  data  published  or  available  shows  us 
that  Van  der  Hoff  could  be  quite  mistaken.  His  paper  on  the  extrapolation  of 
his  function  of  true  stress  versus  strain,  is  based  on  the  "best"  least  squares 
linear  fit  of  measured  stress  at  low  strains.  But  several  points  deserve 
careful  scrutiny  of  all  the  known  stress  versus  strain  data  published  for  low 
values  of  strain  (See  Figs.  25,  26,  27  and  28).  In  particular  we  note  several 
things  from  these  figures. 

1.  That  over  a  large  vange  of  X(from  X  -  1.05  toX  =  1.  60),  the  true 

stress  strain  curve  indicates  a  neo-Hookeian  type  behavior,  i.  e.  the  stress 
strain  curve  is  linear.  We  have  fitted  these  data  to  a  straight  line  over  a  range 

of  X  s  1.05  to  X  :  1.  80.  The  low  variances  and  the  variances  of  the  coerficients 
indicate  a  good  least  squares  fit. 

2.  Often  it  is  more  important  to  indicate  what  a  curve  does  not  do,  rather 
than  what  it  does.  In  all  four  cases  (data  collected  from  Guth,  Shen  and  two  in 
this  study)  the  true  stress  strain  curve  does  not  extrapolate  through  the  origin. 

Even  more  interestingly  the  plots  indicate  that  in  all  four  cases,  there  would  be 
finite  stress  at  zero  strain.  Guth's  data  (Fig.  26)  actually  shows  that  at  very  low 
strains  (1.  01)  there  is  actual  curvature  of  the  true  stress.  This  curvature  is 
drawn  in  dotted  lines  on  the  other  figures,  as  a  trend  rather  than  specific  points. 

It  is  hardly  reasonable  to  suppose  that  all  three  authors  estimated  their  original 

length  with  the  same  sign  of  error.  Rather  if  there  were  errors  in  the  deter¬ 
minations  of  the  rest  lengths,  these  errors  would  be  sometimes  positive  and 
sometimes  negative.  This  is  not  the  case.  AU  samples  even  though  they  are 
different,  show  finite  intercepts  on  the  stress  axis  at  infinitesimal  strains.  In 
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this  context  it  should  be  noted  that  in  3  of  the  analyses  .namely  that  of  Guth  and 
our  two,  the  stress  strain  curves  were  derived  from  force  temperature  data, 
while  Shen's  were  measured  directly.  It  would  be  interesting  to  be  able  to  plot 
the  points  obtained  by  Van  der  Hoff,  before  he  extrapolates  them  to  X  -  1. 

3.  Coming  back  to  Van  der  Hoff's  measurements,  he  postulates  normal  stress- 
strain  behavior,  fits  his  experimental  points  to  the  various  functions,  and  chooses 
the  function  which  shows  the  least  fit  variance  as  the  best  function  which  he  wiil 
extrapolate  to  rest  length.  .This  seems  to  be  a  dangerous  argument  since,  if 
the  stress -strain  behavior  is  different  at  very  low  strains  from  that  at  low  to 
moderate  strains,  it  could  be  incorrect  to  extrapolate  the  rest  length  following 
the  function  describing  the  points  at  moderate  strains.  Guth's  data  (Fig.  26) 
indeed  seem  to  show  chat  the  point  corresponding  to  X  s  1.  01  are  not  on  the 
linear  curve  described  by  the  stress  -  strain  relationship  at  higher  strains. 

The  reason  for  insisting  on  this  seemingly  minor  point  is  that,  as  was  de¬ 
scribed  in  our  Appendix  and  in  Van  der  Hoff's  paper,  the  value  of  the  rest 
length  is  essential  for  deriving  the  behavior  of  the  Mooney-Rivlin  plot,  and 

as  will  be  seen  later,  also  the  fraction  l  it  (inti  *nal  energy  contribution  to  the 
retractive  force),  at  low  values  of  strain;  it  is  there  that  one  would  expect  the 
Gaussian  approximation  to  be  fulfilled. 


4.  A  last  point  to  be  brought  up  in  this  connection  is  the  fact  that  the  de¬ 
parture  from  neo-Hookeian  behavior  is  remarkably  similar  for  all  four  samples. 
Taking,  as  the  measure  for  this  departure,  the  value  of  the  intercept  of  the  true 
stress  at  X  r  1  from  the  least  squares  line,  and  .  p-j  divided  by  the  true 

stress  at  some  arbitrary  value  of  X  within  the  linear  stress-strain  limit  of  all 
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Fig.  29 


Stress-Strain  Curve  -Engineering  -  DAI -22  ,  DAI -23 
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Since  the  rubbers  ere  different  in  the  three  author*'  case*,  tt  it  not  to  be 
expected  thet  the  value*  of  00  /  9lt  50  *hould  be  equivalent,  yet  they  are 
at  least  of  the  same  order  of  magnitude.  Such  is  also  the  ease  for  the  value* 
of  negative  intercept*  on  the  strain  axis,  extrapolated  along  the  least  square  line: 


Guth:  X 

0=  0 


Shan:  X 

9*  0 


this  study: 

DAI -22  X 

arc 


DAI -2  3  x 

<7  ■  0 
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0.  994 


Figure  29  shows  the  engineering  stress  as  a  function  of  strain?  for  both 
our  samples  at  two  different  temperatures.  Using  the  computer  programs  as 
described  previously  in  Appendix  i  ,it  is  possible  to  calculate  the  stress -strain 
behavior  and  any  of  the  derived  plots  for  as  many  temperatures  as  desired. 
Figure  30  shows  the  Mooney -Rivlin  plot  for  sample  DAI -22  and  for  various 
temperatures.  A  detailed  discussion  of  the  constants  2C^  and  2C^  will  be  under¬ 
taken  in  a  later  paragraph. 

Comparing  now  the  Mooney-Rlvlin  plots  of  the  data  taken  from  Guth  (Fig.  31) , 
Shen  (Fig,  32)  and  this  work  (Fig.  33),  we  see  that  in  all  four  cases,  the  Mooney- 
Rivlin  plot  shows  an  upward  departure  from  linearity  as  the  sample  approaches 
its  rest  length.  Even  though  many  precise  determinations  of  the  stress-strain 
curve  at  very  low  elongations  would  be  necessary,  we  do  not  necessarily  believe 
Van  der  Hoff's  assertion  that  at  low  elongations  the  value  of  2C^  vanishes.  How¬ 
ever  we  recognize  our  own  data  taken  around  X  -  1,  is  not  sufficiently  precise 
to  enable  us  to  positively  reject  this  aspect  of  Van  der  Hoff's  work.  Returning 
to  the  Mooney-Rlvlin  plots  of  our  data  (Fig.  30,  33)  (samples 
DAI -22,  DAI -2 3)  the  weighting  function  as  described  in  the  Appendix  I  is  applied 
to  the  weighted  least  squares  linear  fit  for  the  purpose  of  measuring  2C^  and  2C^. 
The  linear  portions  of  the  Mooney-Rlvlin  plots  are  roughly  parallel  as  can  be 
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seen  from  Fig.  34.  The  point*  for  both  DAI  -22  and  DAI -23  are  statistically 
equivalent  and  are  given  for  various  temperatures.  At  25°C,  the  values  of  2C^ 
and  2Cg  are  respectively,  13.  57  •  10®  dynes / cm^  (1.  38  kg/cm^)  and  10.  24  * 

10®  dynes/cm^  (1.  05  kg/cm^).  The  value  of  10^  is  often  less  than  that  of  2C^ 

as  was  mentioned  in  the  Introduction.  In  our  case,  however  there  was  a  finite 
non  sero  value  of  2C^  calculated  for  intermediate  values  of  strain. 

The  temperature  dependence  of  these  coefficients  is  given  below  and 
graphed  in  Figure  35. 

Temp  2Ct  •  2C 

1  2 


• 

c 

10  ®  dynes  W 

10”®  dynes/ 

25 

13.  57 

10.  24 

30 

13.  77 

10.  27 

35 

14.  00 

10.  32 

40 

14.21 

10.  36 

45 

14.  42 

10.41 

50 

14.  63 

10.  46 

55 

14.  83 

10.  51 

60 

15.  04 

10.  57 

TABLE  III  :  Temperature  dependence  of  ZC^  and  2C^ 

ZC,  is  a  linear  function  cf  absolute  temperature.  Comparing  this  table  with 
the  data  obtained  by  Roe  and  Krlgbaum  [  8]  and  T.  L.  Smith  [28],  we  see  that 
2C^  is  appreciably  temperature  dependent,  but  that  2C^  is  practically  temperature 
independent.  For  a  better  view,  the  values  of  2C^  and  2C^  of  Table  III,  normalized 
to  30  C  are  compared  to  the  corresponding  values  of  2C^  and  2C^  at  30°  C  of  Roe 
and  Krigbaum,  and  are  listed  in  the  next  table. 
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T#C 

This 

work 

ZC^ZCi 

at  30°  C 

Roe  and 
Krigbaum 

2C j / 2C j 

at  30°C 

This 

work 

2C2/ 2C2 
at  30°C 

Roe 

Krigbaum 

2C2/2C2 

at  30° C 

30 

1.  000 

1.  000 

1.  000 

1.  000 

40 

1.  035 

1.  032 

1.  010 

1.  009 

50 

1.  065 

1.  062 

1.  019 

1.  028 

60 

1.  094 

1.  093 

1.  028 

1.  042 

TABLE  IV  :  Compare!  son  of  C ^  and  Cg  with  literature  values 

Thus  the  temperature  dependence  of  2Cj  is  the  same  as  that  of  Roe  and 
Krigbaum.  The  temperature  dependence  of  2C  ,  still  small,  seems  to  be 
lar  ,er  in  their  case.  Ours  resembles  more  closely  the  behavior  of  T.  L.  Smith's 
data  [  28l . 

One  of  the  very  important  questions  in  the  theory  of  rubber  elasticity  is 
the  question  of  the  contribution  of  the  internal  energy  to  the  total  elastic  re¬ 
tractive  force.  Up  to  recently  it  was  thought  that  most  of  the  retractive  force 
at  moderate  elongations  came  from  the  decrease  in  entropy  due  to  the  re¬ 
distribution  of ~r*  towards  higher  values,  and  possibly  to  some  chain  alignment. 
However  lately,  it  has  been  shown  that  in  the  case  of  natural  rubber,  something 
like  20%  of  the  force  is  due  to  the  internal  energy  even  at  small  strains.  By 
application  of  Equation  (31)  it  is  possible  to  calculate  fe/f  knowing  (^f/^T) 

P,  -L* 

and  f  vs.  A  at  different  values  of  the  temperature.  Figure  36  shows  the 

values  of  fft/f  for  both  our  samples.  There  is  so  much  scatter  in  the  data  due 

to  the  uncertainty  in  (?tf/&T)  T  ,  that  little  can  be  said  about  the  value  of  f  /f 

p,  j_i  e 

for  low  strains  except  that  there  does  seem  to  be  an  energetic  contribution  to 
the  total  retractive  force  at  low  strains  ,  *»hich  shows  a  decrease  at  higher 
strains,  in  accordance  with  the  literature. 

The  only  available  figures  on  the  dependency  of  f e / f  on  X  at  low  values  of 
strain  come  from  a  paper  to  be  published  by  Shen,  McQuarrie  and  Jackson  [50l. 
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Fig.  38  Internal  Energy  Contribution  from  Shen’s  Data 
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Figure  37  shows  their  values.  While,  for  higher  strains,  all  the  values  in 
recent  publications,  seem  to  agree,  Shan  and  co-workers  at  their  low  strains 
show  a  very  strong  dependence  of  f#/f  on  X .  A  glance  at  Equation  (31) 
shows  us  that  here  again  we  have  a  particularly  sensitive  function  of  X  ,  since 
X  -  1  goes  to  0  as  Xll. 

In  the  following,  we  give  an  analysis  of  the  data  presented  in  tabular  form 
by  Shen  and  co-workers.  Just  as  Van  der  Hoff  showed  the  effect  of  a  small 
change  of  LQ  (the  rest  length)  on  the  Mooney-Rlvlin  plot,  we  can  show  now  the 
effect  of  a  small  change  of  on  the  f#/f  plot.  Knowing  a  given  change  in  Lq 
for  example,  £  0.  01  L.  ,  we  can  recalculate  their  values  of  f  ft,  and  plot  them 
as  in  Figure  38.  This  figure  shows  very  emphatically,  as  did  the  earlier  one 
corresponding  to  the  Mooney-Rlvlin  plot,  the  effect  of  very  small  (1%)  changes 
in  rest  length  on  the  plot  of  fQ/f.  Thus  once  more  we  would  have  to  know  the  rest 
length  with  extreme  accuracy  to  learn  the  true  energetic  contribution  at  very  low 
strains. 

We  can,  now  summarise  this  topic  of  fe/f  by  saying  that,  at  low  to  moderate 
strains,  (X  1.  2  to  2.  5)  the  contribution  to  the  total  retractive  force  of  the  internal 
energy  is  of  the  order  of  15  to  20%,  At  higher  elongations,  (X  >2.  5)  there  is  a 
loss  of  internal  energy  (fe/f  becomes  negative)  due  to  alignment  of  the  chains  that 
are  about  to  crystallise.  At  very  low  elongations  (X<  1-  2)  the  picture  is  not  at  all 
clear.  While,  the  Gaussian  approximation  should  be  best  observed  experimentally, 
according  to  the  data  of  Shen  and  co-workers  the  internal  energy  contribution  rises 
sharply  while  X-*  1. 

Taking  the  data  of  Shen  for  a  Mooney-Rivlin  plot,(2Cj  -4-  2C£)  shows  an 
upturn  at  low  strain  which  Van  der  Hoff  asserts  does  not  exist  if  the  initial  length 
is  very  carefully  determined.  However  one  must  doubt  Van  der  Hoff's  assertion 
of  Hookeian  behavior  (straight  line  extending  from  the  origin  into  the  moderate 
strain  region)  as  discussed  above.  It  remains  then  to  find  a  better  explanation 
of  the  stress  pattern  at  very  low  strains,  especially  the  positive  intercept  on  the 
stress  axis.  Combining  the  evidence  of  Shen's  finding  of  high  internal  energy 
contributions  at  low  strains  with  the  downward  curvature  towards  the  origin  of  the 
stress-strain  plot  we  can  postulate  the  existence  of  a  type  of  superstructure 
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in  the  rubber  that  Imposes  torn*  long  rang*  ordar  at  the  vary  lowest  strain, 
this  auparstructura  holds;  it  must  ba  brokan  up  for  tha  rubbar  to  asauma 
trua  through  going  random  bahavior  at  largar  strains.  This  is  tha  sama  as 
saying  that,  at  vary  low  strains  anargy  must  ba  suppliad  to  distort  tha 
auparstructura.  Onca  tha  latter  has  baan  broken  down,  than  tha  antropic, 
or  cooperativa  rotational,  changas  of  conformation  taka  ovar.  Mora  data 
at  these  low  strains  would  halp  clarify  and  corroborata  tha  pictura. 

A  final  paragraph  may  ba  davotad  to  tha  problem  of  equilibrium  versus 
non -equilibrium  results.  Ever  since  Ciferri  and  Flory  [25]  showed  tha  magnitude 
of  C j  was  influenced  by  tha  time  scale  of  tha  measurements,  much  debate  has 
baan  going  on  as  to  the  nature  of  C Unfortunately,  since  tha  relaxation  of 
force  is  linear  with  log  t,  no  trua  equilibrium  values  have  baan  reported  and 
therefore  tha  interpretation  of  as  an  equilibrium  phenomenon  has  to  be 
postponed.  In  view  of  this,  one  must  also  doubt  that  Van  der  Hoff's  measure¬ 
ments  ware  at  true  equilibrium.  Thus  tha  whole  question  of  C ^  is  still  very 
much  unresolved. 
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IV-  APPENDIX  I 

General  outline  of  programs: 

—  2 1  ■  This  program  indexes  and  computes  the  raw  data  as  it  is 
read  off  the  xv  recorder,  as  force  and  temperature.  There  follows  a 
listing  of  the  program  written  in  FORTRAN  IV,  as  well  as  an  example 
of  the  listed  output.  In  general  each  program  will  have  two  types  of 
output,  one  punched  that  serves  as  part  of  the  input  to  the  next  program 
and  the  other  written  as  a  permanent  record  of  the  calculations. 
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B.  P  391.  This  program  normalizes  the  data,  calculates  the  normalized 
lengths,  calculates  the  various  types  of  stresses  and  fits  the  data  by 
linear  least  squares  analysis.  The  input  to  this  program  is  but  a  slight 
modification  of  the  punched  output  from  P  3  1. 

The  various  types  of  stresses  used  we  defined  as  follows: 

<r  (engineering)  =  F/Area  of  unstressed  sample 

cr (incompressible)  =  F/Area  of  stressed  sample 

where  by  incompressible  we  assume  that  there  is  no  volume  change  upon 
linear  extension.  If  such  is  the  case,  then  where  F  is  the  total  force, 
from  its  geometry,  it  is  easy  to  calculate  the  area  at  any  given  stress  or 
strain. 

AREA  (X  )  =  AREA  (l)/\ 

The  program  also  fits  the  stress  (both  engineering  and  incompressible) 
as  a  function  of  the  temperature.  Both  the  written  and  th«  punched  output 
list  the  least  squares  coefficients,  their  variances  and  covariances.  A 
listing  of  the  subroutines  is  given  in  this  section.  There  is  a  plot  option, 
which  allows  the  checking  of  any  mispunched  data  very  quickly.  The  out¬ 
put  from  P  391  is  the  linear  least  squares  fit  of  the  force -temperature 
data.  Up  to  this  point  no  thermal  expansion  corrections  have 
been  taken  into  consideration. 
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C,  P.  392.  Thin  program  calculate!  the  stress-strain  curvea  at  con- 
atant  temperature  from  the  previoua  force -temperature  curvea  at  con- 
atant  lengths.  With  thia  program  added,  one  now  haa  a  complete  3 - 
dimenaional  picture  of  the  atreaaea  aa  a  function  of  atrain  and  temperature. 

Some  aort  of  deciaion  must  be  reached  as  to  the  order  of  the  fitting 
polynomial.  When  a  testing  program  was  written,  two  important  con¬ 
clusions  became  evident.  Firstly,  aa  the  order  of  the  polynomial  . 
increased,  the  variance  decreased  markedly;  there  ia  a  point  however, 
where,  as  the  order  of  the  fitting  polynomial  keeps  on  increasing,  the 
variance  increases  abruptly.  This  ia  to  be  expected  since,  as  the  fit 
gets  better,  the  residues  get  smaller,  while  raising  the  residues  to  a 
high  power  can  give  number  sufficiently  small  for  the  computer  to  com¬ 
pound  them  with  a  zero  value  (underflow). 


Variance 

2.  14  .  104 
7.35  .  102 
1.97  .  102 
3.22  .  101 
1.06  .  10 1 
3.86  .  103 


Secondly  it  was  found  that  both  the  x's  and  the  y's  to  be  fitted  should 
be  of  the  same  magnitude.  Or,  as  the  magnitude's  change  with  respect 
to  one  another,  the  coefficients  should  change  by  the  same  ratio.  In  the 
case  of  the  higher  order  polynomials,  the  coefficients  were  found  to 
change  by  a  much  larger  order  of  magnitude.  In  this  present  work,  it 
was  decided  to  use  fourth  order  polynomials  to  fit  the  stress-strain  data, 
and  also  to  multiply  the  data  by  a  constant  such  that  x  and  y  were  of  the 


same  order  of  magnitude. 

Once  the  stress-strain  curves  are  stored  in  memory,  one  can 
proceed  to  the  thermal  expansion  corrections  of  stress  and  strain. 

These  corrections  are  small  as  was  mentioned  previously  but  should  still 
be  taken  into  account. 


Finally  the  program  fits  by  least  squares  the  corrected  stress- 
strain  data.  It  is  important  to  note  that  the  weights  given  to  the  points 


fcHhJSAttMtf 


Ill  -  102 


are  the  reciprocal  of  the  variances.  The  propagation  of  variance,  or 
mean  square  error,  is  described  in  DEMING's,  Statistical  Adjustment 
of  Data[48].  In  general  if 


then 

and 

where 


Y  ***+*>  (App.  1:1) 

dy  =  x  da  +  db  (App,  1:2) 

dy^  =  x^  da^  +  db^  +  2x  da  db  (App.  1:3) 

2 

dy  =  variance  of  y 

2  2 

da  ,  db  =  era,  ob  variance  of  a  and  b 
da  db  =  a  ab  covariance  of  a  and  b 


We  are  assuming  no  error  on  x  .  Equation  (App.  I:3)allows  one  to  cal¬ 
culate  the  variance  on  y  (dependent  variable)  knowing  the  elements  of 
the  reciprocal  matrix  and  assuming  the  independent  variable  to  be 
error  free. 
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D,  P  393.  From  the  previous  program,  it  is  now  possible  to  calculate 
the  Mooney- Rivlin  fit  of  the  data.  From [4]  : 

—2—2  =  2Ct  +  2C2  .  £  (App.  I:4) 

The  form  of  this  equation  and  the  importance  of  the  coefficients  2C^  and 
-2C2  are  discussed  in  the  main  text.  However,  it  need  not  be  obvious  that 
any  measurements  at  low  strains  (X.-*  I)  must  be  carried  out  with  extreme 
accuracy.  The  following  sections  will  be  devoted  to  the  sensitivity  of  the 
ordinate  ( <r/  X  )  and  to  the  weight  to  be  given  to  it  any  any  given  value 
of  the  abcissa. 


(App.  1:10) 
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Replacing  the  value  of  d\  ^  from  Equation  (App  IilQ)  into  Equation 
(App  1:7  ),  we  obtain 


(App.  1:11) 


This  may  be  written  as 
dy^  =  Bj  d 0^  + 


(App.  1:12) 


where  B.  and  B  are  functions  only  of  X.  .  There  are  therefore  two  con- 

X  « 

tributions  to  the  total  variance  on  y  ,  the  first  from  the  variance  of  the 
stress  and  the  second  from  the  variance  of  the  unstressed  length.  It 
will  be  seen  in  Appendix  I  that  the  larger  effect  comes  from  the  second 
term.  As  previously,  the  weight  given  to  each  point  in  this  plot  is  equal 
to  the  reciprocal  of  the  total  variance  of  that  point. 

A  careful  scrutiny  of  B  j  and  B2  show  immediately  that  by  far  the 
most  important  factor  is  the  variation  of  the  denominator  as  a  function 
of  X  .  Clearly,  as  X-*  1,  both  denominators  approach  zero,  hence  the 

variance  of  y  approaches  infinite  as  X  approaches  unity.  The  important 
conclusion  of  this  study  is  therefore  that  as  the  elongation  goes  to  unity, 
or  in  other  words  when  the  sample  is  close  to  its  unstressed  length,  the 
uncertainty  of  the  corresponding  Mooney-Rivlin  stresses  tends  towards 
infinity. 
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E,  £  39*.  From  the  punched  output  of  P  392,  it  ia  now  poaaible  to  re¬ 
construct  the  corrected  force -temperature  curves.  It  suffices  to  inter¬ 
polate  the  corrected  stresses  as  is  the  case  of  P  392,  as  a  function  of 
temperature,  keeping  the  strain  a  constant.  A  least  squares  linear  fit 
gives  the  slope  and  intercept  of  the  force -temperature  lines,  which  are 
of  thermodynamic  importance. 
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"FAST -STRETCHING"  of  NATURAL  RUBBER 
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Fig,  39  :  Calibration  of  Linear  Potentiometer 
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III.  EXPERIMENTAL 

For  the  purpose  of  determining  1 1*.  e  temperature  rise  in  the  sample 
upon  elongation,  it  was  necessary  to  embed  a  small  temperature  measur¬ 
ing  device.  This  "thermometer1,  must  have  a  large  thermal  response,  since 

o 

we  are  interested  in  analyzing  temperature  changes  of  the  order  of  0.  01  C. 
Referring  back  to  Part  II,  it  is  now  much  more  obvious  wh^  that  study  had 
been  undertaken.  The  use  of  thermocouple  sensors  requires  very  sensitive 
micro-voltmeters ,  the  usual  sensitivity  of  a  thermocouple  such  as  copper- 
constantan  is  0.  040  mv  per  degrees  Centigrade.  The  measurement  of  0,  001  c; 
would  require  the  measurement  of  4  •  10*^  mv,  which  is  outside  the  realm 

of  normal  type  instrumentation. 

Thermistors  exhibit  several  Interesting  features,  one  of  which  is  a  very 
high  temperature  coefficient  of  resistance.  For  example,  with  the  type  of 
bead  thermistor  that  we  used,  a  sensitivity  of  1  0  mv  per  degree  Centigrade 
is  easily  obtainable.  This  means  a  sensitivity  of  about  250  times  the  magnitude 
that  could  be  obtained  with  a  thermocouple  sensor.  Embedding  the  thermistor 
in  a  -iit  made  in  t  :<s  rubber  with  a  razor  blade,  this  should  have  sittic  c  e<  ; 
upon  the  stress  distr iubution  throughout  the  sample,  even  though  the 
t h er momer e r  head  is  finite,  approximately  0.  03  inches  in  diameter. 

Once  th.e  thermistor  probe  is  embedded,  it  is  held  in  place  by  a  little 
drop  of  rubber  cement  applied  to  each  side  of  the  slit.  By  operating  a  Linear 
Power  Activator  (POLY NOID} ;  m o d e i  m.i m her  0 1  A .  ( Sk i n n e r  Precision  Indus- 
tries.  New  Britain,  Conn.  )  it  is  possible  to  stretch  the  sample  over  a  distance 
of  about  6  inches  in.  less  than  0,  1  secs.  This  is  a  fast  stretch,  ar.d  below 
will  be  referred  to  as  such  (fast -stretch  tests).  As  will  be  seen  in  a  further 
section,  to  label  these  tests  as  being  truly  adiabatic  requires  further  discussion. 

The  increase  in  length  upon  extension  can  be  measured  with  a  Linear 
Motion  Potentiometer .  Such  an  instrument  was  purchased  from  Bourns,  inc. 
Riverside,  California;  model  number  156,  with  a  total  travel  length  of  6.  S  ", 
Knowing  the  rest  length  of  the  sample,  the  stressed  length  from  which  the 
extension  commences,  it  is  possible  to  calibrate  the  linear  potentiometer 
directly  in  cen,  .meters  of  rubber  sample.  This  was  done  in  1  ig.  3b.  It  is 
therefore  easv  to  read  off  the  stressed  length  directly  from  ’•  ig.  3b,  and  !o 
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calculate  the  elongation  of  the  sample  knowing  what  its  unstressed  length 
is.  The  final  stressed  length  is  determined  by  the  position  of  the  rest  against 
which  the  Polynoid  Activator  comes  to.  Small  magnets  hold  the  Activator's 
core  in  place  and  prevent  needless  bouncing  as  the  sample  is  stressed  or 
relaxed.  The  height  of  the  rest  is  controlled  by  a  lab-jack,  which  was 
mounted  in  such  a  manner  that  the  platform  can  move  without  changing  the 
position  of  its  central  torque  axis.  A  hole  was  drilled  through  the  constant 
temperature  chamber,  and  the  stretching  of  the  sample  was  controlled  from 
the  outside. 

All  tests  were  run  at  constant  temperature.  It  was  much  easier  to  analyse 

the  data  if  the  base  lines  were  roughly  horisontal.  A  25  W  white  light  bulb 

connected  to  a  small  powerstat  sufficed  to  keep  the  chamber  at  approximately 
o 

30  C.  A  15  W  fluorescent  bulb  produced  a  very  noisy  Interference  with  the 
thermistor  circuit  and  had  to  be  eliminated.  Figs.  40,  41  show  small  portions 
of  base  lines  the  differences  between  them  is  attributed  to  the  noise  emitted 
by  the  fluorescent  lamp. 

In  the  course  of  these  tests,  a  study  of  the  time  dependency  of  a  stress 
"spike"  was  also  carried  out.  This  was  done  by  activating  the  load  cell  or 
transducer  and  measuring  the  change  in  force  with  time. 

On  the  other  hand,  the  progress  of  relaxation  was  not  followed  at  the 
highest  temperature.  Instead,  the  sample  was  stretched  a  few  times  to  a  ^ 
of  approximately  2.  50,  and  then  hold  there  for  about  12  hours  for  conditioning. 
In  parallel  with  the  measuring  of  the  thermal  effects  of  "fast-stretching, "  we 
will  now  develop  certain  relationships  that  were  implied  in  Part  II.  There  it 
was  shown  that: 
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Part  II:  Equation 
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Thus  while  the  reduced  output,  A«H  ,  at  the  Inflection  point  le  a  function  of 
(a),  the  bridge  ratio:  the  tangent  at  the  inflection  point  It  Independent  of  the  (a) 
ratio.  In  most  cases  to  obtain  any  sensitivity  at  a  given  point,  it  Is  necessary 
to  compensate  or  buca  the  output  st  the  given  point  and  just  measure  the  net 


voltage  of  the  total  circuit.  This  presents  several  Inconveniences,  the  most 
Important  of  which  la  the  difficulty  to  obtain  a  very  stable  compensating  or 
bucking  source.  Looking  at  Equation  11:12,  It  is  obvious  that  Aetl  is  a  function 
of  and  A  which  are  fixed  and  (a),  the  bridge  ratio  ,  which  is  variable.  Thus 
by  changing  the  value  of  (a)  (obtained  by  changing  either  Rj  or  R2  ),  one  can 
supress  the  reduced  output  to  any  value  desired.  Over  a  large  range  of  (a) 

(a  ;  0  to  a  ;  4-  *o  ),  the  value  of  Aelt  varies  from  Tt  4-  0.  5  to  Ti  -0.5, 

r  nr 

and,  as  in  our  case  Ti  c*  0.  1,  varies  approximately  from  0.  6  to  -  0.  4. 

-*r 

There  exists  therefore  a  value  of  (a)  ,  such  that  z  and  for  which  the 

tangent  at  the  Inflection  point  Is  finite  and  Independent  of  (a).  Experimentally 
then,  it  is  possible  to  measure  changes  of  temperature  of  very  small  magnitude 
(0.  001°C). 


While  carrying  out  these  fast -stretch  tests,  it  was  noticed  that  the  thermistor 
output  was  not  returning  to  ltn  original  base  line,  even  after  a  relatively  long 

period  of  time.  From  a  rough  calculation  of  heat  dissipation  through  the  sample, 

It  was  deduced  that  all  of  the  heat  should  have  been  dissipated  within  a  few 
minutes  of  the  stretch  so  that  the  thermistor  could  not  Indicate  residual  stretching 
heat.  From  Fig.  41  it  can  be  seen  that  at  the  end  of  some  3  minutes  (  x  axis: 

50  secs/cm)  the  signals  reached  a  steady  value.  Even  though  the  base  lines  are 
different,  they  are  parallel.  Therefore,  either  the  thermistor  is  responding  to 
some  other  variable  than  temperature,  or  there  is  some  process  taking  place  in 
the  sample  which  is  exothermic,  it  is  also  very  important  that  the  applied 
potential  across  the  Wheatstone  bridge  be  a  constant.  A  small  drift  in  E  can 
give  a  considerable  change  in  base  lines. 

e 

We  examined  the  latter  hypotheses  first,  but  since  neither  relaxation  nor 

crystallisation  could  be  seriously  considered  (stress  relaxation  work  convinced  { 

us  that  under  these  conditions  the  relaxation  was  minimal,  while  at  small  I 

1 

elongations  of  the  order  of  100%  there  could  not  be  any  crystallisation  phenomenon  j 
accounting  for  an  exothermic  process),  we  turned  our  attention  to  the  possible  I 
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TRUE  STRESS  AS  CALCULATED  FROM  ASSUMING  THE 
MATERIAL  INCOMPRESSIBLE 
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response  of  the  tharmiator  to  aoma  othar  varlabla.  Plotting  tha  ravaralbla 
difference  in  baaa  Unaa  aa  a  function  of  atraln  (Kig.  42)  ,  wa  found  that, 
naglacting  aoma  acattaring,  a  straight  lina  could  ba  drawn  through  the  point 
at  which  wa  atartad  to  atrateh.  Thla  obaarvatlon  combinad  with  tha  fact  that, 

If  atratchlng  took  placa  slowly,  tha  diffaranca  In  baaa  lines  observed  was  the 
same  as  for  a  fast  atrateh,  lad  us  to  conclude  that  tha  tharmiator  was  pressure 
sanaitlva  and  that  tha  plot  of  (Fig,  42)  amounts  to  a  measure  of  true  ati  ass  aa 
a  function  of  strain,  A  careful  comparison  of  the  true  stress  aa  a  function  of 
strain,  (Fig,  43)  showa  a  linear  relationship  between  thla  and  tha  course  of 
tha  tharmiator  base  line.  Thua  monitoring  tha  fast  stretching  experiment 
by  a  thermistor  leada  to  tha  combined  maaauramant  of  tha  atraaa  and  tha 
temperature  rise.  Unfortunately,  there  is  too  much  scatter  In  the  points 
(Fig.  42)  to  enable  us  to  use  thla  technique  to  measure  tha  atraaa  in  tha  rubber 
directly.  However  due  to  tha  differences  in  time  response  of  pressure  and 
temperature  on  tha  tharmiator  (temperature  time  response  la  practically 
Immediate),  the  temperature  rise  upon  stretching  la  given  by  the  total  peak  , 
while  the  incremental  preseure  le  given  by  the  differences  In  base  lines. 


A.  Calibration  of  Thermistor  4 


Using  a  LEEDS-NORTHRUP  CO.  ,  Wheatstone  bridge,  catalogue  number 
4760,  in  conjunction  with  a  large  2  volt  WILLARD  battery  and  a  Dc  miliivoltmeter 
from  SENSITIVE  RESEARCH  Instrument  Corp.  ,  Mount  Vernon,  New  York, 
we  measured  the  resistance  R  of  the  thermistor  at  a  series  of  temperatures. 
Treating  the  R-T  data  as  shown  in  Part  II,  it  was  readily  possible  to  calibrate 
the  thermistor  (Appendix  II). 
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Run  # 

recorder 

rubber  (cmi) 

final 

MC 

"Ap"cn 

1 

30.  6 

11.  8 

1.  62 

0.  028 

i' 

0.  017 

2 

0.  028 

2* 

0.  014 

3 

40.  2 

13.  6 

1.  87 

0.  073 

0.  55 

3' 

0.  064 

0.  60 

4 

0.  065 

0.  60 

4* 

0.  058 

0.  60 

5 

50.  0 

15.  5 

2.  12 

0.  105 

0.  90 

5' 

0.  068 

1.  00 

6 

_ 

0.  105 

0.  95 

6' 

0.  078 

1.  00 

7 

62.  2 

17.  9 

2.45 

0.  170 

l.  30 

7* 

0.  160 

1.  30 

8 

0.  162 

1.  30 

8' 

0.  157 

1.40 

9 

77.  8 

21.  0 

2.  88 

0.  232 

1.  65 

9' 

0.  224 

1.  70 

10 

0.  226 

1.  55 

10* 

0.  207 

1.  40 

11 

68.  0 

19.  1 

2.  60 

0.  loS 

1.  00 

11' 

0.  168 

1.  20 

12 

0.  170 

1.  10 

12* 

0.  165 

1.  25 

13 

54.  5 

16.  4 

2.  25 

0.  129 

0.  85 

13* 

0.  109 

0.  85 

14 

0.  126 

0.  85 

14' 

0.  120 

0.  90 
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TABLE  V  Cont'd 


L  of 


Run  # 

r scorder 

rubber  (cms) 

X 

final 

A  TCr 

"A  p" 

15 

43.  6 

14.  3 

1.  96 

0.  081 

0.  60 

15' 

0.  078 

0.  60 

16 

0.  086 

0.  60 

16' 

0.  07  3 

0.  60 

17 

36.  0 

12.  •: 

1.  72 

0.  061 

0.  50 

17' 

0.  056 

0.45 

18 

0.  067 

0.  60 

18* 

0.  056 

2.  50 

emt 
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RESULTS 


A  typical  sample  of  natural  rubber  cut  out  from  the  same  sheet  as  those 
used  in  other  experiments  (number  DAI -26)  was  employed  In  the  experiment 
described  in  detail  below.  It  had  the  following  dimensions: 

Zero  length  =  80  mm,  width  -  4.02  mm,  thickness  -  l.  90  mm,  length  ; 

at  rest  -  73  mm,  area  -  7.638  *  10“^  cm^ 

The  operating  temperature  was  2S.  25  0  C,  and  the  thermistor  circuit  set  at 
(Fig.  44): 

E  =  l  volt,  :  1430a,  R2  -  999.9A  R4  :  1412  a,  a:  1000 

The  stretching  of  the  sample  starts  from  a  slightly  elongated  position  (X  -  1,  10) 
so  as  to  minimize  errors  attending  the  "rest"  state  and  also  to  reduce  the  effect 
of  the  relaxation  spike.  The  whole  problem  will  be  dwelt  upon  in  detail  below  in 
the  Discussion  section  but,  briefly,  what  seems  to  happen  is  that  when  the  rubber 
is  rapidly  brought  back  by  our  sample  holder  from  the  stretched  to  the  unstretched 
condition,  its  own  retraction  speed  is  slower  than  the  movement  of  the  solenoid, 
and  its  shape  becomes  momentarily  buckled  during  the  relaxation  to  aero  force. 

The  sample  returns  then  to  its  original  shape  and  elongation  of  1. 1  slowly  and 
thus  returns  to  a  state  of  slight  sti  is.  This  return  is  better  defined  than  one  to 
complete  rest. 

Fig.  39  shows  the  linear  calibration  used  for  determining  the  increase  in 
length  of  the  sample  upon  stretching.  Table  V  shows  the  final  elongations  to  which 
the  rubber  was  stretched,  as  well  as  the  differences  in  base  lines  and  changes 
in  temperature. 


Figure  45  shows  the  increase  and  decrease  of  temperature  as  a  function  of 
elongation,  remembering  that  the  original  elonption  from  which  rubber  U 
stretched  is  1.  10. 
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ig,  47  :  Relaxation  of  spikes  at  54°  C  using  bucking  potential 
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The  experiments  of  fast  -  stretching  approximate  adiabatic  conditions.  If 
the  stretching  of  the  sample  Is  performed  sufficiently  fast,  then  the  heat  developed 
on  account  of  the  loss  of  entropy  {  dQ  rev"  Td.S)  measurable  as  the  temp¬ 
erature  rise  of  the  sample.  Retracting,  the  sample  turns  colder,  as  the  entropy 
returns  towards  its  maximum  value  at  rest.  In  the  procese  of  retraction,  in  the 
Ideally  adiabatic  experiment,  one  should  be  able  to  measure  a  hyeterleie  in  the  rate 
of  return  of  the  sample  to  the  unstressed  state  since  A  S  is  positive,  the  reaction 
endothermic  and  yet  no  heat  available  from  the  environment.  In  the  case  of 
atomic  gases,  which  would  expand  in  analogy,  the  increase  in  entropy  would  be 
Immediately  balanced  by  corresponding  cooling.  In  the  case  of  rubber  the 
additional  number  of  conformations  can  ba  assumed  only  by  drawing  on  the  energy 
reservoir  of  rotational  and  vlbra  ionaal  dagraea  of  freedom  which  rnesni  that 
while  the  rubber  becomes  more  random  it  also  becomes  leas  mobile.  This  ia 
the  same  process  which  causes  thegts  atoms  to  cool  but  cooling  in  the  case  of 
rubber  means  lower  mobility  and  thereby  a  slower  return  to  its  rest  length.  In 
other  words  ,  for  rubber  to  return  to  its  original  shape  with  the  same  speed  with 
which  it  could  be  extended  would  require  an  influx  of  heat  or  a  not  fully  adiabatic 
experiment. 

It  is  then  Interesting  to  notice  the  time  effect  on  the  force  (Fig.  46)  when  the 
sample  is  being  stretched,  while  all  three  variables,  namely  the  rise  in  temp¬ 
erature  of  the  sample,  the  length  to  which  the  sample  is  being  stretched,  and  the 
total  force  applied  on  the  sample,  art  being  plotted  simultaneously.  There  seems 
to  be  a  relaxation  of  the  force  to  a  reasonably  steady  value  after  several  minutes 
after  the  sample  has  been  stretched  fast  (1/10  second)  to  a  constant  length.  This 
phenomenon  was  much  more  clearly  established  when  a  counter  e.  m.  f.  (bucking  ' 

force)  of  approximately  12  mv  was  used  to  measure  the  total  e.  m.  f.  of  13  mv.  ' 

(Figs.  47,  48).  In  order  to  elucidate  the  nature  of  these  force  peaks,  two  functions 

were  examined.  Firstly,  the  effect  of  the  pre -elongation  from  which  tha  sample  1 

is  stretched,  and  secondly  the  Influence  of  temperature.  From  Figures  49 
and  50  it  is  clear  that  for  stretches  to  the  same  final  length,  the  lower  the  initial 
length,  the  steeper  is  the  relaxation  from  the  maximum  height  of  the  epi  ke. 

Figure  49  shows  the  relaxation  curves  for  the  elongating  samples.  The  equilibrium  j 
force  attained  is  the  same  whatever  the  initial  elongation,  as  it  should  be.  One 
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can  calculate  a  relaxation  iime  for  the  recovery  from  the  overshoot  force  which 
turn*  out  to  be  the  ahorter  the  more  the  sample  has  been  pre-stretched.  This  dif¬ 
ference  in  stretching  distance  arises  of  course  from  the  fact  that  different  elong¬ 
ations  were  ob  'ned.  We  find  for  the  sample  stretched  from  rest  length  to  a 
X  _  2.  00,  a  relaxation  time  of  about  150  minutes  for  the  slow  relaxation  range 
of  the  sample  stretched  from  X  -  1.  50  to  2.  00,  As  might  be  expected  from  the 

sharpness  of  the  spike  marking  the  overshooting  of  force,  there  is  also  a  very 
fast  relaxation  regime  which  for  all  the  runs  is  of  the  order  of  r  r  1.5  minutes, 
interestingly  enough,  the  relaxation  process  lor  this  internal  local  over  extension 
of  the  sample  has  only  two  relaxation  steps.  On  retraction  the  samples  start 
from  the  same  strain,  arrive  at  the  different  elongations  so  that  equilibrium  forces 
at  the  end  of  the  experiment  are  correspondingly  different.  Again  an  overshooting 
is  observed,  now  towards  the  low  side  of  the  force.  In  fact  the  sample  can  be 
observed  temporarily  to  buckle  so  that  it  is  obvious  that  it  does  not  retract  with 
its  spontaneous  speed  but  is  being  pushed.  The  recovery  to  the  equilibrium  force 

again  differs  for  the  3  samples  but  the  relaxation  times  are  difficult  to  establish 
because  the  time  when  the  sample  assumes  its  equilibrium  length  are  hard  to 
ascertain. 

To  explain  these  observations  it  is  postulated  that  the  overshooting  of  the 
force  at  the  end  of  the  stretching  is  due  to  the  high  speeds  of  the  test,  too  high 
for  the  chains  to  follow  the  motions  in  an  affine  manner  and  the  net  effoct  is  an 
increase  in  labile  entanglements,  that  act  as  temporary  additional  cross-links 
and  raise  the  modu.us.  When  the  motion  stops,  the  force  decreases  steadily 
because  of  relaxation  and  disappearance  of  the  network  of  labile  crosslinks. 

Stating  it  differently,  the  chains  are  impeded  in  their  extension  and  aligning  by 
preexisting  entanglements  and  internal  viscosity.  Thus  one  can  understand  that 
the  relaxation  of  the  force  in  the  experiments  producing  the  greatest  elongation 
should  have  the  fastest  decay  since  in  this  case  more  labile  crosslinks  (shorter  - 
T  )  have  been  created.  Pre-stretching  the  sample  to  an  appreciable  elongation 
means  that  the  chains  are  more  equilibrated  towards  the  final  strain  and  there  ie 
leas  additional  entanglements  on  further  stretching. 

During  the  retraction  process  complications  of  molecular  motion  should  be 
less  due  to  entanglement  processes  and  more  to  the  fact  that  the  internal  friction 
may  not  allow  the  chains  to  return  to  their  random  conformations  during  the 
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Fig.  52:  Eff«ct  of  Temperature  upon  retraction  of  rubber 
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instrumental  time  interval.  This  lag  will  be  enhanced  by  the  largely 
adiabatic  character  of  the  retraction  which  permits  the  chains  only  to  a 
small  extent  to  acquire  the  necessary  Influx  of  enthalpy  from  the  surround¬ 
ing  s  and  has  to  rely  largely  on  conversions  from  the  internal  energy  reservoir. 
The  picture  is  one  of  a  substantial  number  of  molecules  not  recovering  their 
rest  length  fast  enough,  the  sample  stays  longer  than  if  it  would  contract  in 
less  than  the  instrument  time  and  in  fact  it  can  be  seen  to  buckl  a  and  bow 
during  the  fast  push  back.  In  this  context,  it  is  interesting  to  follow  the 
effect  of  temperature  on  these  phenomena.  Figures  51  and  52  show  the 
action  of  temperature  on  the  relaxation  of  the  force  as  a  function  of  time. 

Figure  51  demonstrates  that  at  higher  temperature  the  force  relaxes  faster. 
This  is  to  be  expected,  especially  on  the  basis  of  our  picture  of  labile  entangle¬ 
ments.  One  can  see  also  on  Figure  51  that  the  results  are  completely  re¬ 
producible,  (runs  1  and  2)  and  that  the  points  lie  on  the  same  relaxation  curve. 
Figure  52  shows  the  comparable  behavior  for  the  retraction  relaxation  curves. 
It  would  be  interesting  to  examine  the  phenomenon  of  time  dependent  moduli 
as  affected  by  the  speed  of  stretching  more  quantitatively.  Ordinarily  the 
relaxation  times  of  a  well  vulcanised  rubber  are  too  short  to  become  noticeable 
during  ordinary  stretching  experiments.  In  our  case  we  have  a  rather  weakly 
cross-linked  rubber  which  permits  entanglements  to  become  more  effective 
and  to  form  time  dependent  cross-links  which,  as  the  experiments  show, 
exhibit  much  longer  relaxation  times. 

Finally  the  specific  heat  of  our  rubber  at  constant  length  and  pressure 
was  calculated  using  Equation  63  : 


and  integrating  this  expression  between  L  and  L,  assuming  the  approximation  I 

®  t 

that  C ^  p  is  independent  of  T,  which  is  acceptable  because  of  the  small  AT  » 
we  obtain  that 

c  :  t  j(i f/»r)p,L  dL 

pV  AT 
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Fig.  51  Engineering  Stress  -  Temperature  Coefficient  as  a 
function  of  Uicorrected  Normalized  Length. 
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Thla  expression  can  be  further  simplified,  by  setting 

X=  L/L0 

or  dX  =  dL/L 


and  since  V  -  V  -  L  A  ,  then 
o  o  o 


where  A  :  initial  area 
o 


with  f 


T  IX  * 

S  =  Ai 

L.  p  _ 

p  AT 

as  the  stress  per  unit  undeformed  area. 


Calculation  of  the  integral  of  the  force -temperature  coefficient  la  done 
by  the  use  of  a  planimeter.  Figure  53  shows  the  plot  of  (Af  /d  T) 

\  p  p  L 

as  a  function  of  X,  All  elongations  start  at  X  z  1.  1  for  reasons 
that  we  have  already  discussed. 


X  final 

area  of  (tf/dT)  L 
cm2 

(*/<m 

10-3  P'L 

P 

eat. 

1.  50 

0.  031 

1.  225 

0.91 

1. 75 

0.  072 

2.  840 

0.91 

2.  00 

0.  121 

4.  780 

0.  92 

2.25 

0.  179 

7.  060 

92 

2.  50 

0.  247 

9.750 

0.93 

2.  80 

0.  351 

13.  880 

0.93 

TABLE  VI:  Values  of  (%/VT) 

P» 

L  mnd  P 

estimated 

Application  of  £quatior/69)  permits  one  to  calculate  C  as  mentioned  before. 

O  P 

The  specific  heat  at  300  k  is  given  below 
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X  final 

\J  eng  p,  *- 

1  -5 

*  10 

A  T 

C  ,  ca’ 

1.  50 

3.  675 

0.019 

—  P,  -L.  — — 

0.  505 

1.75 

8.  520 

0.  045 

0.495 

2.  00 

14.340 

0.  080 

0.465 

2.25 

21.  180 

0.  120 

0.455 

2.  50 

29.  250 

0.  165 

0.452 

2.  80 

41.  640 

0.  225 

0.  470 

TABLE  VII  :  Values  of  C  ^  at  different  elongations 


Thus  the  specific  heats  calculated  from  the  adiabatic  temperature  rise  and 

from  the  force  temperature  measurements  are  quite  similar  to  those  of 

natural  rubber  measured  by  calorimetric  means  [  51  ]  1.  e.  have  a  value  of 
o  o 

0.  480  cal/g  C  at  300  K.  The  agreement  Is  deemed  good  considering  the 
fact  that  there  is  appreciable  scattering  of  the  measurement  of  ^Tatnabatic 
and  also  the  fact  that  we  were  engaged  here  in  an  indirect  method  of 
measurement. 

The  specific  heat  of  the  rubber  as  expected,  decreases  with  increasing 
elongation.  Our  numbers  should  be  taken  as  approximate  rather  than  quantita¬ 
tively.  The  decrease  in  specific  heat  observed  is  10%  .  Bekkedahl  and 
Matheson  ^5l1  have  shown  that  there  is  a  decrease  of  some  8%  in  the  specific 
heat  of  partly  crystalline  rubber  versus  the  amorphous  sample.  When  chains 
stretch  and  align,  there  is  a  loss  of  some  degrees  of  freedom.  In  the  random 
conformation,  the  molecules  have  3  vibrational  degrees  of  freedom,  and  some 
of  the  possible  3  rotational  .  When  they  align  and  or  start  to  crystallise,  there 
Is  a  freezing  of  vibrational  degrees  of  freedom  along  the  long  dimensions  of  the 
chains  and  also  of  some  rotational  motions  expressing  Itself  as  a  net  decrease 
in  specific  heat.  It  would  be  very  interesting  to  measure  the  specific  heat  of 
various  elastomers  at  different  elongations  and  temperatures.  These  results 
should  then  be  compared  to  the  results  obtained  from  direct  calorimetric 
measurements,  and  x-ray,  I.  R,  and  birefringence  data. 
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Run 

T°C 

R  JL 

0 

T  C 

R-ft 

T°C 

RJI 

1 

3  3.  99 

1740 

29.  96 

2018 

25.  91 

23.  48 

2 

33.  94 

1744 

29.  92 

2023 

25.  90 

2350 

3 

33.  91 

1744 

29.  89 

2025 

25.  88 

2353 

Ave. 

33.  95 

1  ’42 

29.  91 

2022 

25.  89 

2350 

TABLE  VIII:  Resistance  as  a  function  of  temperature 

B. 

Determination  of  8  30.  00°  C 

T°C 

t°k 

1  /T  1/°K 

R  -A. 

lo«10  R 

$  °  K 

33.  95 

307.  10 

0.  0032563 

1742 

0.  2410 

3446 

29.  91 

303.  06 

0. 0032996 

2022 

0.  3058 

3385 

25.  89 

299.  04 

0. C033440 

2350 

0. 3211 

TABLE  IX:  Determination  of  9  30.  OO^C 

where 

*  :  2. 

3026  (logjQ  Rj  - 

lo«l0 

V 

i 

1/Tj  -  1/T0 

and  at 

29.  91°C, 

3415°  K, 

R  - 

2022  JV 

i 
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Calibration  of  Thermistor 
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D.  Checking  of  Therm  stor _ £  4 


^eu  ;  T1/Rl  -  0.500  for  a  :  1000 
A e 1 1  ;  -  0.4113  mv/mv 

Using  E  :  100  mv,  Rj  :  100A  R2  :  0.  1 A  end  R4  :  1412  A 


Aen  -  -  0.4137  mv 

11  *  lire 


*•'“  :  «Ti*  -  ».2 

t,2 


-  9.  002  '  10"3  mv  /°C 


30.  33 
30.  44 
30.  56 
30.  7  1 


See  Figure  54  ,  which  show*  the  calibration 
of  the  thermistor  actually  imbedded  in  the  sample! 
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DILATOMETRY  OF  NATURAL  RUBBER 
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EXPERIMENTAL 


The  change  In  total  volume  of  a  typical  rubber  during  stretching  as 
we  have  seen  previously,  is  very  small,  of  the  order  of  10“^  cc/cc  [36,  37 

38]  .  Our  basic  design  of  the  dilatometer  and  the  method  of  testing  resembles 
closely  the  methods  used  by  Dr.  Buckely  at  the  Esso  Research  Laboratories. 
Figures  55,  5c.  5  <  show  sketches  of  the  aluminum  dilatometer  machined 
by  R.  Parla  at  the  Brooklyn  Polytechnic  Institute.  Following  Buckely's 
claims,  we  planned  first  to  place  the  dilatometer  in  an  air  bath  controlled 
to  4  0.  01°  C;  it.  was  thought  then  to  be  easy  to  obtain  a  control  of  It  0.  001°C 
in  the  dilatometer  cavity  on  account  of  the  large  heat  sink  constituted  by  the 
mass  of  metal.  Constancy  to  0.  001  C  was  thought  necessary  to  ensure  that 
thermometer  effects  would  be  an  order  of  magnitude  smaller  than  the  volume 
changes  due  to  the  stretching  of  the  samples.  Buckely  and  his  group 
used  Nichrome  wire,  commercially  available,  and  showed  that  the  variation 
in  thickness  of  the  wire  was  negligible,  so  that  the  construction  of  the  di¬ 
latometer  could  be  based  on  a  design  with  pull  from  the  outside  .  As  will  be 
shown  in  the  Results  section,  this  was  proved  to  be  correct.  Unfortunately 
other  essential  features  of  the  Esso  dilatometer  proved  to  be  unworkable,  and 
our  dilatometric  technique  had  to  be  extensively  changed.  In  particular,  Buckely 
and  co-workers'  use  of  water  as  the  confining  liquid  like  others  quoted  before, 
again  supported  our  original  idea  of  working  at  3.  97°C,  at  which  the  co¬ 
efficient  of  thermal  expansion  of  water  is  nil,  had  to  be  abandoned,  since  water 
was  found  to  sweii  rubber,  even  over  short  periods  of  time  [52]  . 

It  was  also  noticed  that  our  natural  rubber  became  discolored  upon  long 
Immersions  in  water  and  a  detailed  analysis  of  the  behavior  of  this  elastomer 
in  various  liquids  was  undertaken.  Fig.  58  shows  the  percent  weight  change 
as  a  function  of  time  of  natural  rubber  in  various  liquids.  T  52"1  While  no 
claim  is  made  as  to  the  change  in  volume  upon  immersion,  these  tests  do  show 
a  distinct  and  reproducible  change  in  weight.  A  truly  acceptable  non-inter¬ 
acting,  confining  fluid  for  the  dilatometric  studies,  must  not  cause  any  measur¬ 
able  uptake  over  a  substantial  period  of  time.  Various  salts  were  added  to 
water.  Their  effect  on  the  "swelling"  is  shown  in  Fig.  58.  Many  other  fluids 
were  tried  such  as  fluorocarbons,  other  known  non-solvents,  mono,  di-,  and 
trifunctional  and  alkali  salts.  Glycerine,  pentanediol  and  Li  salts  were  the 
only  non-swelling  agents  found  for  natural  rubber.  It  was  rationalised  that  Li 
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ross-sectionai  view  of  constant  Temperature  bath  tor 


■alts  war*  ln*rt  because  as  a  raault  of  thalr  large  aolublllty  and  becauae  of 
th*  high  charge  to  volume  denalty  of  the  Li-salts  they  hold  their  eolvating 
water  moat  tightly.  Li-chlorides  aa  well  ae  nitrate*  aeemed  to  have  no 
ewelUng  cffecte.  Rather  than  pentanediol  or  glycerine,  which  were  too 
highly  viacoue  and  hygroacoplc  under  our  condition*  it  wee  decided  to 
uae  the  concentrated  nitrate  aolutlona  becauae  of  better  metal  paaeivation. 

For  a  period  of  over  200  hour*,  there  wae  no  detectable  weight  pickup,  but 
unfortunately  theae  aalta  had  a  corroding  effect  on  the  aluminium  block. 

Varioua  teata  were  run,  and  it  waa  ahown  that  a  3  to  5  mil  depoalt  of  Nickel 
metal  plated  electrolyticaily  on  the  aluminium  aeemed  to  protect  the  latter 
from  corroeion.  All  holder*,  ecrewa  and  other  acceaaorlea  were  alao  nickel 
plated  or  made  of  etalnles*  ateel.  A  nechrome  wire  waa  uaed  to  elongate 
the  eample.  The  aample  h cider  waa  of  the  type  (aelf-tightenlng  wedge*  with 
rubber  pade)  aa  that  uaed  in  the  force-temperature  meaaurementa  (Fig  13). 

A  detail  of  the  pad*  la  alao  shown.  They  were  bonded  to  the  wedges  with  an 
epoxy  resin  and,  as  the  sample  ia  stressed,  are  kept  in  place  by  tightening 
of  the  wedges.  Again,  this  type  holder  was  found  to  be  very  efficient  at 
moderate  to  high  elongations. 

The  design  of  the  dllatometer  ia  featured  further;  2  rubber  gaskets  form 
a  tight  seal  with  the  glass  plate  (thickness  of  5  mm)  to  which  a  female  joint 
is  attached,  A  specially  constructed  doubly  inverted  teflon  seal  was  designed 
to  allow  passage  of  the  nichrome  wire  through  the  block,  without  leakage. 

Fig.  57  shows  a  detailed  croaaectlonal  view  of  this  special  gasket.  The 
interchangeable  capillaries  attached  to  a  fitting  male  joint  can  be  inserted  in 
the  dllatometer  fitting.  It  waa  found  that  the  thermometer  effect  of  the  LiNO-j 
solution  in  the  cavity  required  0.  0001°C  to  ensure  base  'line  stability.  A 
constant  temperature  bath  controlling  to  better  than  0.  001°C  was  then  designed 
(Fig*  59,  60). 

To  this  end,  a  few  basic  principles  concerning  thermal  control  had  to  be 

observed  or  elaborated.  Firstly  all  individual  factors  must  be  isolated  and 

tested  separately.  A  thermistor  bead  was  used  to  ascertain  the  fluctuations 

in  temperature  in  various  areas  of  the  bath.  Careful  checking  with  Beckman 

thermometer  in  conjunction  with  the  thermistor  tested  for  random  departures 

>4  o 

from  the  base  line  of  no  mors  than  about  5*10  C.  It  was 
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essential  to  hold  the  temperature  within  0.  001°C  for  several  days.  To  this 
end,  a  lightening  heavy  duty  stirrer  was  placed  in  each  corner  of  the  bath. 

The  thermistor  control  probe  was  hooked  up  to  a  Bayley  proportional  con¬ 
troller,  model  123  ,  Bayley  Instrument  Company,  Danville,  California  to 

sense  and  anticipate  the  changes  in  temperature.  This  unit  is  guaranteed 
£ 

to  tO.  001  C.  In  view  of  the  good  insulation  a  glo-quartr  immersion  100 
W  heater  was  the  only  heat  source  for  the  overall  bath. 

The  evaporation  of  the  water  from  the  bath  was  the  single  biggest  dis¬ 
turbance.  It  was  found  necessary  to  enclose  the  water  by  a  surrounding  air 
bath  and  to  heat  the  air  about  3°C  above  the  water  temperature  which  was 
kept  at  34° C  and  to  control  the  air  temperature  by  a  thermlstemp  on-off 
contoller,  model  171.  Slabs  of  polystyrene  foam  provided  thermal  lneulatlon 
from  the  room,  and  a  squirrel-cage  fan  was  used  to  mix  the  air.  Since  the 
lightening  stirrers  were  running  quite  hot,  the  air  temperature  was  maintained 

Q 

at  approximately  37  C  by  blowing  in  cool  air  with  a  fan.  To  maintain  a  water 
temperature  over  a  period  of  days  of  1  0.  0005° C  greatest  care  must  be  taken 
in  handling  the  whole  system.  Two  holes  were  cut  out  of  the  front  polystyrene 
slab  and  rubber  gloves  inserted.  Further  to  perform  the  stretching  of  the 
rubber,  two  axles  were  guided  through  the  walls  of  the  bath  gasketed  with  0 
ring  seals.  By  turning  the  shafts,  the  samples  can  be  stretched  without 
disturbance  of  the  thermal  equilibrium. 

The  solutions  of  60  g  of  LiNO^  /100  cc  of  H  O  were  refluxed  for  about 
24  hours  and  filtered  through  mllipore  filters.  The  rubber  samples  were 
placed  into  the  dilatometer,  all  gaskets  tightened  with  a  torque  wrench,  and 
the  concentrated  LiN03  solution  poured  into  the  dilatometer  cavities.  The 
latter  were  thendsgasnd  thoroughly  by  applying  a  vacuum:  since  small  bubbles 
could  cause  variations  in  capilllary  height,  due  to  either  atsmospheric  and 
or  thermal  changes,  orders  of  magnitude  larger  than  the  effect  we  were 
looking  for.  For  example,  a  change  in  pressure  of  0.  01  atm.  on  a  bubble 
of  .  1  cc,  would  cause  the  level  in  the  capillary  to  change  by  1  cms  while  the 
expected  change  in  volume  by  etretching  the  rubber  is  of  the  order  of  1  mm  ! 
Small  thermal  changes  create  of  course,  similar  changes  in  the  height  of  the 
capillaries. 

After  the  dilatometer  is  placed  into  the  bath  gross  thermal  control  is 
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obtained  by  using  the  thermietemp  on-off  controller.  After  a  period  of 
approximately  12  hours  to  ensure  thermal  equilibrium,  the  excess  LiNO^ 
solution  is  wiped  off  the  top  of  the  capillaries,  and  the  bath  Is  allowed  to 
cool  a  few  hundreths  of  a  degree.  The  proportional  controller  Is  then  switched 

on,  and  the  whole  thermal  system  is  allowed  anew  to  come  to  equilibrium. 

The  height  In  the  capillaries  is  plotted  as  a  function  of  time  to  establish  the 
base  lines.  Generally  after  a  horizontal  or  a  slowly  varying  base  line  has 
been  reached,  both  samples  ware  stretched  simultaneously  and  the  changes 
in  volume  plotted  as  a  function  of  time. 

In  view  of  the  paramount  Importance  of  maintaining  constant  temperature 
of  at  least  I  0.  001°  C  over  long  periods  of  time  as  stressed  above  the  true 
temperatures  of  the  bath  and  in  the  dllatometrlc  cavities  were  also  recorded. 

To  this  end,  thermistor  probe  was  placed  near  the  Beckman  in  the  bath,  and 
the  readings  from  the  Beckman  as  well  as  the  changes  in  height  of  the  capillaries 

Q 

followed.  The  Beckman  thermometer  showed  a  stability  at  34  C  of  at  least 
0.  00 1*  C  over  a  span  of  several  hours,  while  the  thermistor  probe  and  the 
height  of  the  liquid  in  the  capillaries  indicate  a  stability  better  than  5  •  10"^°C. 

Testing  for  the  degree  of  leak  proofedness  of  the  dllatometer  wires, 
pulling  them  with  no  sample  present  showed  no  movement  in  the  capillaries. 

We  can  thus  assume  also  further  that  there  are  no  variatlonss  in  the  diameter 
of  the  wires. 

The  height  of  the  liquid  $  oluanna,  in  the  capillaries 
(precision  bore)  was  calibrated  as  a  function  of  degrees  Centigrade  by  plotting 
the  height  in  the  capillaries  versus  temperature  read  on  the  Beckman.  We 
ascertained  an  average  value  of  59.  5  cm/degree  Centigrade  for  the  capillaries 
of  diameter  of  0.  0383  cm,  and  a  volume  calibration  of  1.  15  •  10-3  cc/cm. 

The  samples  used  were  natural  rubber  cut  from  the  same  sheet  as  used  in 


the  other  tests. 

length  at 

Thickness 

Width 

Area 

Sample 

rest  (mm) 

(mm) 

10Z  cmZ 

DAI -24 

49 

1.  90 

3.  96 

7.  524 

DAI -25 

49 

7.  83 

4.  11 

7.  521 

m  -  us 


After  the  stability  of  the  ay  stems  was  ascertained  both  samples  (DAI >24 
and  DAI -25)  ware  atretched  by  the  aame  amount  at  the  eame  time.  Aa  both 
aamplea  were  from  the  aame  aheet,  one  would  have  expected  the  aame  results 
from  both  cavitiea.  Unfortunately  ,  before  a  aufftetent  no.  cauUbt  tested,  we 

ran  into  aome  alow  leakage  problema  of  the  gaaketa  and  were  able  to  get 
acceptable  data  for  juat  2  runa.  The  trouble  ia  almoat  certainly  due  to  corroalon 
of  the  metal  of  the  gaaket  grooves.  The  deciaion  to  place  the  dllatometer  in 
water  inatead  of  in  an  air  bath,  and  the  rubber  into  aalt  solution  inatead  of  into 
pure  water  had  to  be  taken  long  after  the  dllatometer  construction  ao  that 
aluminium  turned  out  to  be  the  wrong  metal  to  use.  For  further  work  the 
dllatometer  should  be  conatructed  of  a  stainless  steel  resistant  to  salt  corrosion, 

and  the  flat  gaskets  be  made  out  of  teflon  to  minimise  their  swelling. 
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(MINS)  Fig.  62  :  Dilation  run  No. 
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Fig,  64  ;  Relative  Dilation  for  sample  DAI -25 
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n  RESULTS 

In  view  of  the  limitations  discussed  above,  we  present  the  diiatometric 
data  more  to  show  trend  and  magnitude  rather  than  as  a  quantitative  result, 
figures  61  and  62  show  the  behavior  of  the  capillaries  as  a  function  of  the 
elongation  of  the  samples  and  the  time  for  which  they  were  held  at  a  given 
elongation.  As  can  be  seen ,  the  base  line  is  very  stable  for  the  rear  cavity 
(DAI -25)  and  slightly  increasing  for  the  front  cavity  (DAI -24).  Thus  it  would 
be  difficult  to  subtract  the  ill -defined  base  line  to  obtain  the  dilation  effect  on 
sample  DAI  -24.  f  igures  63  and  64  show  the  same  results  expressed  in 
changes  of  volume  (AV/Vq  •  10"^)  as  a  function  of  normalized  length. 

There  can  be  no  doubt  that  the  volume  of  the  rubber  increases  from  the 
vary  first  stretching,  but  that  the  amount  of  increase  is  very  small. 


Ill  -  169 


Calibration  Wheel 

Sample:  DAI  -24 

Length 

(mm) 

-  0.  2394  cms/tooth  VQ  -  0. 

Cap.  -  ] 

X  ^V(cm)  AV(cc) 

3700 

l.  15  •  10'3 

^v/v 

cc/cm 

cc/cc 

49 

1 

0 

0 

0 

49  -A  23.  94  = 

72.  94 

1.49 

0.  015 

1.72  ■ 

M 

o 

1 

yj\ 

4.  65 

•  io-5 

49  4-  47.  88  s 

96.  88 

1.  98 

0.  045 

5.  20  • 

10-5 

1.40 

b* 

O 

1 

49  A-  71.82  = 

120.  82 

2.  46 

0.  110 

1.  26 

O 

l 

* 

3.40 

o 

1 

-4 

.  .  _-4 

49  -A  95.76  ; 

144. 26 

2.  96 

0.  120 

1.  38  ’ 

10 

3.  74 

10 

49  A-  105.  33  - 

154.  33 

3.  24 

0.  140 

1.  56  • 

■ 

o 

4.  40 

•  io-4 

TABLE  X: 

Ditatometry,  First  run 

on  sample  DAI -24 

Calibration  Wheel 

=  0.  2394  cms/ tooth 

V  - 
o  - 

0.  3700 

Sample:  DAI  -24 

Cap.  - 

I 

o 

l A 
^4 

f-4 

% 

3  cc/cm 

Length  (mm) 

X 

^V(cm) 

^V(cc) 

AV/V^cc 

/cc 

49 

1.  00 

0 

0 

0 

53.  8 

1.  10 

0 

0 

0 

58.  6 

1.  19 

0.  005 

5.70  ’  10"6 

1.54  • 

IO'5 

63.4 

1.  29 

0.  010 

1.  15  •  10’5 

3.  12  • 

IO-5 

68.  2 

1.  39 

0.  015 

1.72  *  IO"5 

4.65  - 

• 

10“5 
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0.  015 
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*-* 
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I 

U1 
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o 

1 

vn 
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7.  75  1 
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1.  69 

0.  045 
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o 
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Length  (mm) 

X 

AV(cm) 

AV(cc) 

AV/** 

cc/ cc 

87.  34 

1.  79 

0.  065 

7.45  • 

IO-5 

2.02 

•  io"4 

91.  14 

1.  89 

0.  075 

8.  60  ‘ 

io-5 

2.  32 

■  IO*4 

96.  88 

1.  98 

0.  080 

9.  20  • 

io'5 

2.48 

■  IO-4 

TABLE 

XI: 

Dllatometry,  Second  run  on  sample  DAI 

-24 

Length  (mm) 

X 

A  V(cm) 

A  V(cc) 

AV/V 

cc/ cc 

72.  94 

1.  49 

0.  030 

3.45  ' 

IO"5 

9.  30 

•  10-5 

96.  88 

1.  98 

0.  060 

6.  90  • 

10-5 

1.  89 

•  10‘4 

120.  82 

2.46 

0.  060 

6.90  ’ 

10"5 

1.  89 

*  IO"4 

144.  76 

2.  96 

0.  030 

3.45  • 

io-5 

9.  30 

•  10‘5 

154.  33 

3.  24 

-  0.  040 

•  4.  60  • 

IO"5 

-  1.  22 

•  10*4 

TABLE 

XII: 

Diiatometry,  First 

run  on  sampl 

e  DAI- 

25 

Length  (mm) 

X 

AV(cm) 

AV(cc) 

.  *V/X 

cc/cc 

49 

l.  00 

0 

0 

0 

53.  8 

1.  10 

0.  010 

1.  15  ■ 

-5 

10 

3.  1  • 

IO"5 

58.  6 

1.  19 

0.  020 

2.  30  * 

IO-5 

6.  2  • 
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63.  4 

1.  29 
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8.8  • 
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68.  2 

1.  39 
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3.45  * 
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8.  8  • 
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72.  9 

1.49 

0.  030 

3.45  * 

IO’5 

8.8  • 

IO'5 

77.  7 

1.  59 

0.  030 

3.45  • 

10‘5 

8.  8  ' 

IO-5 
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Length  (mm) 

X 

A  V(cm) 

AV  (cc) 

AV/Vn 

cc/cc 

82.  5 

1.  69 

0.  030 

3.45  * 

10"5 

8.  8  • 

10-5 

87.  3 

1.  79 

0.  045 

4.  80  * 

10“5 

1.  30  * 

IO-4 

91.  I 

1.  89 

0.  040 

4.25  * 

I0-5 

l.  15  • 

i 

o 

96.  8 

1.  98 

0.  040 

4.25  * 

io-5 

1.  15  • 

i 

o 

TABLE  XIII:  Dilatometry,  Second  run  on  sample  DAI -25 
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m.  DISCUSSION 

One  of  the  most  Important  aspects  of  this  study  of  dilation  was  the 
swelling  of  natural  rubber  In  different  liquids.  Water  itself  is  deflnately 
a  swelling  agent,  up  to  3%  weight  change  is  recorded  after  an  immersion 
of  several  days,  li  this  weight  uptake  were  translated  entirely  Into  & 
volume  change,  this  would  correspond  to  a  AV/Vq  of  3  •  10  which  is  two 
orders  of  magnitude  above  the  expect  :  change  from  dilation  due  to  the 
Isotropic  tension. 

Thus  one  of  the  first  considerations  is  to  question  any  results  in  the 
literature  that  give  the  dilation  of  natural  rubber  from  experiments  carried 
out  in  water,  almost  the  only  kind  of  experiment  used  to  date.  Holt  and 
McPherson  [161  and  Wolstenholme  T35]both  measured  the  dilation  in  water, 
although  the  former  do  mention  that  for  longer  periods  of  immersion  they 
used  mercury.  However  we  found  that  water  uptake  was  certainly  not 
negligible  at  even  small  immersion  times.  For  the  resuits  obtained  by 
hydrostatic  weightings  T36,  381,  one  should  ask  oneself  how  much  water  is 
taken  up  immediately  by  capillary  action  into  the  rubber  voids.  A  fraction 
of  a  percent  would  be  still  an  order  of  magnitude  above  the  expected  mechan¬ 
ical  dilation.  According  to  Grove's  resuits  [53l;  water  is  not  automatically 
eliminated  as  dilatometrlc,  fluid  because  much  of  the  imbibed  liquid  does  not 
seem  to  change  the  overall  volume.  But  even  if  the  overall  volume  did  not 
change  during  weight  pick-up  of  liquid,  the  liquid  f mbihed  must  respond  to 
the  stretching  in  several  ways: 

1.  by  flowing  out  of  the  larger  voids  and  pores,  it  will  make  the 
dilatometrlc  volume  of  the  rubber  upon  stretching  seem  larger. 

2.  parallel  to  this  effect,  the  liquid  that  has  gone  into  the  molecular 
interstices  will  only  reappear  in  the  bulk  of  the  liquid  very  slowly, 
and  hence  restricts  the  motion  of  the  rubber  chains  and  prevent 
accomodation  to  the  new  length.  This  will  tend  to  decrease  the  Poisson 
ratio  and  make  a  contribution  to  the  apparent  large  volume  Increase  upon 
stretching. 

Lowering  the  activity  coefficient  of  the  water  by  addition  of  highly  soluble 
salts,  we  found,  as  mentioned,  gave  negligible  swelling  after  several  days  of 
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testing. ,  Since  K  or  Na  salts  were  not  so  effective,  one  has  to  deduce  that 
the  Li  with  its  stronger  capacity  to  bind  water  molecules  quite  strongly, 
inhibits  water  swelling  of  the  rubber  network  at  concentrations  where  there 

is  no  free  water  left  r 541  ,  This  in  itself  was  a  very  interesting  result. 

Before  examining  and  discussing  the  experimental  data,  It  might  be 
of  Interest  to  analyze  on  the  molecular  level  why  an  isotropic  body  expands  upon 
simple  tension.  In  the  case  of  liquids,  which  relative  to  the  applied  stresses 
are  incompressible,  all  forces  immediately  cause  motion.  The  packing  in  a 
liquid  is  loose  enough  (the  order  is  short  range)  so  that  the  molecules  need 
not  penetrate  over  a  geometrical  barrier,  as  there  are  enough  voids  to  take 
up  any  motion  and  stress  relaxation  is  practically  instantaneous.  In  the  case 
of  an  elastic  solid  or  a  glass,  the  molecules  require  more  than  fluctuation 
space  to  move  past  surrounding  denser  clusters.  During  recoverable 
oeformation  the  molecules  need  not  move  from  one  trough  to  another,  but 
merely  move  up  the  energy  quadient  and  then  fallback  into  the  old  position  on 
removal  of  the  stress.  If  the  molecules  really  moved,  we  would  no  longer  be 
in  the  elastic  region,  but  irreversible  plastic  flow  would  take  place.  Hence, 
as  the  molecules  tend  to  "rise"  on  their  neighbors,  they  Increase  their  inter- 
molecular  distances  and  as  a  result  there  is  a  net  increase  in  volume.  I*  urther, 
as  the  molecule  is  in  a  potential  energy  well  increasing  its  distance  to  some 
of  its  neighbors  must  increase  also  its  potential  energy.  But  potential  energy 
is  free  energy,  and  hence  whenever  there  is  a  change  in  stored  energy  due  to 
distortions!  stress,  there  must  be  a  dilation  in  volume.  An  interesting 
corollary  is  that,  when  in  the  case  of  a  swollen  network  the  molecules  can 
move  by  the  way  of  the  voids  present  in  the  swelling  liquid  there  can  be  an 
elastic  storage  upon  stretching  for  which  one  would  not  expect  any  change  in 
volume.  In  general ,  though,  one  can  say  that  any  isotropic  elastic  body  will 
undergo  dilation  due  to  the  negative  isotropic  stress  component  in  tensile 
deformation. 

The  four  runs  that  were  carried  out  (see  Tables  X,  XI,  XII  and  XIII 
also  Figures  61  and  62)  even  if  not  quantitatively  »*«»Hable,  show  unequivocally 
that  natural  rubber  exhibits  indeed  volume  dilation  upon  small  elongations. 

At  the  higher  elongations,  the  supposedly  identical  samples  begin  to  differ. 

We  are  not  in  a  position  at  this  time  to  explain  the  differences,  and  more  work 


would  have  to  bo  carriad  out  with  a  now  dllatomator.  A  pootlblo  explan¬ 
ation  might  bo  that  tho  motion  of  tho  wires  enhance  oamotic  leakage 
or  promotoa  gaa  nueloatiai  following  corrosion.  As  for  tho  observod  volume 
changes,  wo  conclude  that  tho  volume  expansion  la  very  small,  corresponding 
to  a  Poisson  ratio  of  0.499.  That  Is,  rubber  does  Indeed  behave  at  first 
almost  like  a  true  liquid,  and  volume  expansion  occurs  l  cally  and  near  cross 
lime  points.  In  our  weaxly  cross- linked  rubber  this  must  be  expected  to  be 
a  small  disturbance.  More  tightly  vulcanised  or  partly  filled  rubbers  aught 
to  exhibit  more  dilation  even  initially. 

Eventually,  once  the  experimental  technique  Is  sufficiently  reliable,  It 
is  planned  to  measure,  amongst  others,  the  volume  effects  of  stress  re¬ 
laxation,  the  dliational  reversibility  of  the  stress-strain  cycle,  the  volume 
dilation  at  low  to  moderate  strains  to  establish  the  nature  of  the  Poisson 
coefficient  of  the  rubbers,  and  finally  the  kinetics  and  density  changes  at 
higher  elongations  due  to  crystallisation  .  It  will  be  particularly  worthwhile 
to  study  these  same  parameters  In  the  case  of  non-crystalllslng  rubbers. 
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IV.  APPENDIX  III 

A.  Change  in  volume  upon  uniaxial  tension  -  Theory 

Consider  a  cylindrical  isotropic  material,  of  radius  R0  and  of  length 
Lq.  The  initial  volume  is  given  by: 

VQ  :  tfR2o  L0  (App.  11:1) 

at  a  given  length  L, 

V  :  ffR2  L  (App.  11:2) 


where  R  and  L  are  the  radius  and  the  length  respectively  at  any  given  time. 


Av 

0 

> 

i 

> 

it 

(App.  II:  3) 

where 

V 

=  "(R0  -A  AR)*  (Lc  -A  aL) 

(App.  11:4) 

and 

Vo 

-  trR2  L2 
o  ° 

thus 

A  V 

=  !>(R0  -A  AR)2  (Lq  A-  AL) 

-  t  R2  L  J  In  R2  L 

n  °  O  ° 

(App.  11:5) 


l  Poisson  ratio  defined  in  X  strains 

We  will  define  the  X  Poisson  Ratio,  as  being: 

-  AR/  *0 

AL/Lq 


V  = 


(App.  11:6) 
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« 
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and  if  furthermore 

C  «  AL/Lq  (App.  11:7) 

than  AR/Ro  .  -  (V  (App.  11:8) 

App.  II,  Equation  (5)  can  be  written  as: 

AV/V o  ;  (I  -A  AR/Rq  )2  (1  -A  AL/Lq)  .  ! 

=  (1  A  2  Ar/ro  -A  (AR/Rq)2  )  (1  -A  AL/Lq  )  -  I 

:  1  -A  AL/L  -a  2  AR/R  •  AL/L  -A  (AR/R  )2  -A 

u  °  o 

(AR/Ro)2  Al/Lq  -  1  -A.  2  AR/R0 

:  (*2cv*  2tt''  <  A  e2^  -A  e3  (App.  11:9) 

App.  II,  Equation  (9)  can  be  written  in  two  forms: 

AV/V  -  C  (1  -  2  -A  C2  (»^  -  2  -Ac3!/2  (App.  11:10) 

o 

A  V/VQ  -  (Al/(-2f2  -  2  0  A^I^At3)  (App.  II:  10') 


These  equations  (10,  10')  are  exact  expressions  giving  the  change  of  volume 
of  an  Isotropic  sample  deformed  under  simple  tension. 

a.  Case  of  an  incompressible  body 

If  Av  =  0,  then  the  body  by  definition  is  incompressible.  Setting  AV 
equal  to  zero  in  App.  II:  Equation  (10'), 


vZt  C  C  2  (c  -A  1)1-  2  v[  c  «  -A  1)1  -A  c  :  0  (App.  11:11) 

solving  for  V  : 

vi  •=  <  -a  1)  (c  2  (c  -A  I)2  -  c2  (c  -A  1)  •  c  )1/Z 


t2«A  1) 

z  c  (c  -A  1)  c  (c  2  -A  2c  -A 


c2 


1/2 

0 


f  2  (c  -A  1) 
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<App.  II:  14) 
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What  arc  the  valuta  o£  tha  tangents  to  this  function? 

-3/2  -1/2  , 

vi  -i  T  l  (c  -A  1)  1  -  rl  -  (r  -A  1)  If 


-3/2 


*  f  (f  -A  1)  "  *-  ^  [1  -  (c  -A  irl/21 


(App.  11:15) 


Applying  Hospital's  rula 
At  ‘  ;  0 

(V'i)  : 

C  -  0 


-3/2  -5/2  -3/2 

«  -A  1)  -A  c<«3)  («  -A  1)  -  (c  -A  1) 

2 


2  C 


-5/2  -5/2  -7/2 


;  -3  (f  -A  l)  -3 

2  2 

(C  -A  1) 

-A  15  c(c  -A  l) 

4 

-5/2 

3  (C  -A  1) 

2 

2 

2 

:  -3  -3  -A  3 

2  2  2 

2 

<5 

II 

o 

1  1 

•►1  U) 

(App.  11:16) 

at  f  ;  00  >(v 't) 

:  0 

(App.  11:17) 

f  - 


Thus  one  conclusion  is  immediate,  and  that  is  that  at  very  small  strains, 
the  value  of  the  incompressible  Poisson  ratio  approaches  the  predicted  value 
of  0.  5,  but  along  a  non  aero  tangent.  At  small  strains,  the  value  of  Poisson' a 
ratio  is  not  0.  5  for  incompressible  bodies. 
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i.,,  Ci>»  of  a  compressible  body 
Rearranging  App,  II:  Equation  (10'), 

v2r<2(f-A  1)1-  2V[<(c-A  1)  1  -A  f  -  AV/V0  :  0  (App.  11:18) 


Ir  nnalogy  with  the  compressible  body,  we  shall  choose  the  negative  root. 

,  2  , 

v  z  c  (  c  -A  1)  -  (  c  2  (  f  -A  1)  -  (  1  (c  -A  1)  <  c  -  A  V/VQ)  ) 

<  2  <  c  -A  1) 


1/2 


v  i  (  c  -A  1)  -  (  (  c  -A  1)  (t  -A  1  -  c  -A  AV/Vq  )  ) 


c  (  «  -A  l) 


ill 

Vs  1  [1  -  (V/V0  (  c  -A  1)  ]  (App.  11:19) 

e 


1  rom  App.  II:  Equation  (19)  it  would  be  possible  to  solve  for  one  of  the  thre< 
variables  knowing  the  two  others.  To  enable  one  to  solve  this  equation  easily 
we  will  use  the  reduced  variables. 


Let  us  call  V/V0  -  d  (dilation) 

(  C  -A  1)  -  d  ;  S 

1/2 


then 


*v-  1  “<  A _ ) 

d  -t  S 


1  -  .  1 


1/2 


1  -A  S/d 


) 


Tf 


then 


x  ;  S/d 

y  ;  <  V 

y  -  l  -  (i  -A  x  ) 


-1/2 


(App.  11:20) 

(App.  11:21) 
(App.  11:22) 
(App.  11:2  3) 
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A  two-dimensional  plot  ofy  ;  f(x)  can  be  drawn  up,  and  the  Poisson  ratio 
of  a  compressible  body  as  a  function  of  the  change  In  volume  (d)  and  a  function 
of  the  elongation  (e)  can  be  calculated. 

Z.  Poisson  ratio  defined  In  natural  strains 


If  the  Poisson  ratio,  -  In  ( 1  -A  AR/R0) 


In  (1  -A  AL/Lq) 


(App.  11:24) 


v  :  ff  (R„  -A  AR)  (L„  -A  AL) 


Lo 


r  (1  -A  AR  )  (1  ~b  Al_  ) 


(App.  11:25) 


and  if 


(1  -A  AR  )  -  (  1  -b  AL  ) 


(App.  11:24) 


1-2  V 


(1  -b  AL  ) 


(App.  11:26) 


AV  -  V  -  1  :  (  1  -b  «) 

V*  *  Vo 
o  ° 


1  -  2  v 


(App.  11:27) 


Contrary  to  the  expression  for  AV/V^  obtained  by  using  the  Poisson  ratio 
as  defined  in  X  strains,  this  expression  for  AV/Vq  is  always  zero  for  a 
Poisson  ratio  of  0.  5. 


Thus  an  isotropic  incompressible  body  has  a  Poisson  ratio  as  defined  by 
natural  strains,  always  equal  to  0.  5  and  this  value  of  V  ^is  independent  of 
the  elongation  c  • 
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B.  FI  of  of  Poisson  ratio*  a>  a  function  of  dilation  and  elongation 


From  App.  11:  Equation  (19)  it  it  possible  to  calculate  the  Poisson  ratio 
for  any  given  dilation  and  any  given  elongation.  This  equation  has  been  set 
up  on  the  computer,  and  Figures  (65,  66,  67)  show  the  Poisson  ratio  for 
various  ranges  of  elongation  and  dilation  . 

1  igure  (67)  shows  the  figures  for  both  large  elongations  and  large  dilations 
Here  it  Is  quite  apparent  that  in  the  case  of  Poisson  ratio  defined  from  X  strains, 
a  material  of  such  a  constant  equal  to  0.  5,  would  show  a  strong  contraction 
in  volume  upon  elongation.  Also  at  <  z  200%,  one  would  actually  have  no 
volume  left  at  all  I  Figure  (66)  shows  intermediate  Poisson  ratio  as  defined 
by  X  strains  as  a  function  of  small  values  of  elongation  and  dilation.  There 
it  can  be  clearly  seen  that  the  tangent  at  the  origin  (c  ;  0)  for  V  ^  -  0.  5  is 
not  equal  to  zero.  Thus  even  at  infinitesimal  strains,  for  a  Poisson  ratio  as 
defined  here  equal  to  0.  5,  the  dilation  (actually  a  contraction)  is  finite. 

Figures  (68,  69)  show  the  Poisson  ratio,  as  defined  for  various  natural 
strains,  as  a  function  of  elongation  and  dilation.  In  this  case,  negative  values 
of  dilation  (or  contraction  of  volume)  are  not  possible  for  <  0.  5.  Also 
if  the  material  Is  truly  incompressible,  then  the  Poisson  ratio  is  0.  5  for  all 
elongations. 
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CONCLUSION 

J  Of  the  three  different  experiment*  carried  out  on  Natural.  Aubber,  the 
force  versus  temperature  measurement  have  verified  that  the  decay  of  stress 
in  stretched  elastomers  is  approximately  linear  with  logarithm  of  time  for 
a  considerable  period,,  as  was  found  by  Kraus  and  Moeaugembx  fM]  r'&ent 
[27]  ,  and.  Smith  [W]  .  A  viscoelastic  mechanism  for  this  decay  according  to 
a  spectrum  of  Maxwell  terms,  can  be  ruled  out  as  the  basis  for  this  form  of 
relaxation  because  of  the  nature  of  the  lime  dependence.  As  a  further  difficulty 
small  local  defects  leading  to  crack  propagation  by  stress  concentration  and 
crystallisation  of  the  sample  complicate  this  relatively  simple  but  unexplained 
form  of  relaxation.  Thus,  there  is  a  natural  difficulty  in  the  way  of  de¬ 
termining  when  the  criterion  of  equilibrium  has  been  adequately  fulfilled.  In 
our  experiments,  the  sample  was  considered  fully  relaxed  when  the  Increment 
of  relaxation  over  the  total  span  of  testing  was  less  than  one  percent.  '  Other  ^ 
authors  simply  indicate  without  further  elucidation  that  they  waited  for  the 
decay  of  stresses  to  ha  "negligible". 

Over  the  ranges  of  temperature  (30  C  -  60°  C)  and  elongations  X  -  (1.  00 
-  3.00)  studied,  the  force  temperature  curves  were  reproducible,  re¬ 
versible  and  linear.  A  series  of  computer  programs  were  written  in  order 
to  correct  the  raw  data  for  the  thermal  expansion  of  the  rubber.  (Before 
correction,  the  slope  of  the  force  temperature  data  is  negative  at  low  values 
of  strain,  and  becomes  positive  at  higher  values;  this  phenomena  is  known  as 
the  thermo-elastic  inversion,  first  pointed  out  by  Meyer  and  Ferri  [55]  and 
corrected  for  by  Anthony,  Caston  and  Guth  [46]  .  This  thermo-elastic 
inversion  simply  shows  that  at  the  lower  elongations  the  thermal  expansion 
of  the  rubber  by  increasing  its  length  at  constant  stress,  reduces  the  tension 
at  constant  length.  In  other  words  at  very  low  strains  the  reduction  of  tension 
by  thermal  expansion  exceeds  the  increase  of  tension  with  temperature  to  be 
expected  from  the  kinetic  theory  of  elasticity.  The  thermo-elastic  inversion 
point  is  the  strain  at  which  these  two  effects  balance  exactly.  The  thermo¬ 
elastic  inversion  influences  also  of  course,  the  adiabatic  temperature  change 
in  elastomers,  which  if  (3f/3T)_  ,  remains  uncorrected  will  be  calculated  to 
have  a  negative  sign  at  very  low  strains  (Fig.  23).  Unfortunately,  our  ex¬ 
periments  were  not  sufficiently  accurate  at  these  low  strains  to  observe  this 
inversion. 


The  Mooney-Rivlin  parametric  equation  a*  was  stressed  in  the  Intro¬ 
duction  has  no  molecular  basis,  and  really  describes  a  purely  hypothetical 
body.  It  was  also  shown  that,  even  though  several  authors  [22,  50]  have 
done  extensive  work  on  the  dependency  of  both,  of  the  Mooney -Rivlin 
equation  and  of  fQ/f  (the  internal  energy  contribution  to  the  total  retractive 
force  )  upon  the  elongations  at  low  values  of  strain,  Roe  and  Krigbaum  [18] 
feel  that  precise  data  gathered  from  the  low  strain  values  would  cast  addi¬ 
tional  light  upon  the  validity  of  the  Gaussian  approximation  ,  which  should  be 
valid  as  long  as  the  chains  are  only  slightly  stretched.  Unfortunately,  the 
behavior  of  the  samples  at  these  low  strains  might  well  suffer  from  a  lack 
of  well  defined  response  to  a  stress.  At  moderate  elongations  the  network 
is  deformed  as  a  whole,  and  we  are  truly  seeing  average  properties  taken 
over  a  large  population,  but  at  very  small  strains,  many  fewer  chains,  or 
small  sections  of  the  network,  may  participate  in  the  deformation,  and  the 
averages  become  accidental  functions  of  minor  disturbances. 

Moreover  as  shown  in  the  Appendices,  it  is  mathematically  obvious  that 
at  very  low  strains  the  errors  of  the  measurements,  be  they  on  the  Mooney- 
Rivlin  strain  (  X  -  X  or  the  internal  energy  contribution  (  X  ^  -  1)"^  ,  must 
become  very  large  i.  e.  a  small  uncertainty  on  X  is  greatly  magnified  as  X 
approaches  unity.  The  available  data  at  low  strains,  particularly  that  ob¬ 
tained  from  Anthony,  Caston  and  Guth  [46]  show  in  all  cases  that,  if  the 
actual  stress  is  plotted  as  a  function  of  strain,  a  straight  line  is  obtained  that 
does  not  pass  through  the  origin.  This  we  have  attributed  to  a  type  of  super¬ 
structure,  that  first  must  be  broken  up  (requiring  energy  to  do  so),  before 
the  material  exhibits  truly  rubber  elastic  response.  This  fact  would  tend  to 
negate  Van  der  Hoff's  [22l  ascertion  that  for  very  low  values  of  strain,  the 

term  of  the  Mooney-Rivlin  equation  is  nil.  Van  der  Hoff  defines  his  rest 
length  (critical  because  of  the  sensitivity  of  both  and  fe/f  to  its  value  at 
low  strains  )  by  assuming  Hookeian  behavior  and  placing  the  origin  into  the 
straight  extension  of  the  line  of  his  points.  We  do  not  dispute  a  linear  type  of 
response  at  moderate  strains  (from  X  s  1.2  to  1.8)  but  we  do  question  that  this 
linearity  can  be  presumed  to  extend  into  the  region  of  still  lower  strains. 

With  the  aid  of  a  weighting  factor  it  was  possible  to  calculate  the  values 
of  the  Mooney-Rivlin  parameters  (2C^  and  2C£)  for  the  combined  data  of  the 
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two  samples.  Cj  depends  linearly  on  T  as  predicted  by  the  simple  kinetic 
theory,  while  was  found  to  be  relatively  insensitive  to  changes  in  T. 

Again,  it  is  difficult  to  read  the  data  pertaining  to  C^.  In  our  analysis  ,  all 
data  at  values  of  X  smaller  than  1.  25  were  discarded  as  being  error  prone. 

It  (its  similarity  be  deduced  that  f e / f  must  depend  strongly  upon  the 
rest  length  at  low  elongations  and  that,  in  fact,  Shen's  data  tend  to  show  a 
rapidly  increasing  value  of  C2  in  agreement  with  others  in  the  literature 
and  in  contrast  to  Van  der  Hoff's  data.  Small  changes  in  rest  length  indeed 
changes  the  aspect  of  fe/f.  drastically  from  an  increasing  contribution  of 
internal  energy  at  the  low  strains,  to  a  decreasing  one,  depending  on  the 
chosen  value  of  the  rest  length.  Cne  might  sum  up  the  outcome  of  this  study 
with  respect  to  and  C2  into  the  following  way  with  reference  to  the  idealized 
diagram  where  the  stress  is  directly  proportional  to  the  Mooney-Rivlin  strain 
(X  -  1/X2)  ,  Referring  to  thiB  linear  extension  graph  the  Gaussian  line  rep¬ 
resents  the  final  approximation  with  all  of  its  well  known  oversimplifying 
assumptions.  The  Mooney-Rivlin  curve,  cutting  off  after  the  second  term,  is 
a  good  second  order  approximation.  According  to  this  investigation  C 2  is 
not  zero,  probably  not  even  in  the  limit  of  X-»  1.  Thus  linear  back  extrapolation 
to  the  origin  are  not  justified,  though  an  experimental  verification  is  extremely 
difficult  because  of  the  indetermination  of  the  magnitude  of  deformation  (zero 
length)  at  small  elongations  and  possibly  also  because  of  the  fluctuations  of 
structures  in  undeformed  rubbers  themselves,  and  because  as  soon  as  these 
fluctuations  became  less  important  on  account  of  larger  deformations,  the 
rubber  has  become  anisotropic  and  ordinary  continuum  theories  can  no  longer 
apply.  To  judge,  by  the  thermal  inversion  point,  this  may  become  true  at 
much  smaller  elongations  than  hitherto  thought.  The  first  deviations  from 
Gaussian  behavior,  as  encompassed  by  the  Mooney-Rivlin  theory  and  indirectly 
by  experiment,  are,  similar  to  many  other  cases,  probably  best  understood 
in  terms  of  Van  der  Waals  type  attraction  or  repulsion  forces,  and  here  also 
as  anisotropic  volume  terms.  There  may  be  also,  but  not  decisively  con¬ 
tributions  from  non-equilibrium  states.  While  none  of  these  suggestions  is 
entirely  new,  it  is  believed  that  as  matter  of  their  emphasis,  confirmation  of 
certain  previous  information,  and  interelation  of  these  and  previous  findings, 

the  present  is  a  new  aspect  of  rubberlike  behavior  at  low  and  moderate 
elongations. 
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The  result*  obtained  from  dilatometry  are  encouraging.  Preliminary 
tests  show  that  the  samples  do  dilate  upon  elongation,  and  that  the  dilation  is 
of  the  expected  magnitude  .  of  10  cc/cc.  The  Poisson  ratio,  defined  in  terms 
of  natural  strains,  turns  out  to  be  approximately  0.  499  which  is  very  ciose  to 
showing  that  our  rubber  is  incompressible,  experimentally  the  difficulties 
so  far  could  not  yet  be  overcome.  The  water  bath  in  which  the  di'atometer  is 
immersed  must  be  controlled  to  better  than  0.  001°C  over  relatively  long 
periods  of  time.  Secondly  the  dilatometer  must  be  leakproof  better  than  one 
part  in  a  million  and  stay  quite  uncorroded.  Thirdly,  the  swelling  of  the 
rubber  upon  immersion  of  the  sample  must  be  less  than  1  part  in  1000  and 
change  little  with  temperature  or  state  of  strain.  We  could  fulfill  the  first, 
and  practically  also  the  third  condition,  but  failed  so  far  in  the  second. 

Before  investigating  the  thermal  effects  of  fast  stretching,  it  was  necessary 
to  derive  a  linearized  output  of  the  Thermistor  as  the  fourth  leg  of  a  Wheat¬ 
stone  Bridge  as  a  function  of  temperature.  The  reduced  open-circuit  output 
potential  was  found  to  be  an  inflected  function  of  temperature;  knowing  the 
resistance  and  the  temperature  coefficient  of  the  resistance  at  a  specified 
point  allowed  one  to  chooBe  the  position  of  the  Thermistor  Function  at  inflexion. 
The  slope  of  the  tangent  as  well  as  the  value  of  the  inflection  point  were 
accessible  to  theroretical  calculations. 

The  use  of  a  thermistor  bead  as  a  temperature  measuring  device  presents 
many  advantages  (its  size,  its  sensitivity,  and  its  fast  response  to  thermal 
changes).  Unfortunately,  in  this  particular  case  it  is  also  pressure  sensitive 
and  this  effect  must  be  corrected  for  in  the  total  response.  Since  the  pressure 
response  appears  to  be  much  slower  than  the  thermal  response,  the  original 
peak  may  be  taken  to  be  wholly  due  to  thermal  changes  in  the  rubber.  With 
the  help  of  the  temperature  changes  during  elongation  and  retraction  thus 
measured  and  of  the  previously  determined  force-temperature  curves,  it  was 
possible  to  calculate  the  specific  heats  of  the  stretched  elastomer,  independ¬ 
ent  of  temperature  and  at  various  stages  of  strain.  It  was  found  that  the 
specific  heat  is  a  decreasing  function  with  increasing  elongation.  A  result  not 
unexpected  in  view  of  the  effective  aligning  of  chains  and  therefore  loss  of 
entropy.  Compared  with  the  literature  these  results  are  found  to  agree  quite 

closely  with  the  few  available  other  data. 
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Stress -strain  data  from  isothermal  stepwise  :itretching  experiments  are 
usually  difficult  to  evaluate  because  it  is  so  hard  to  ascertain  whether  one 
deals  truly  with  equilibrium  data.  Therefore,  not  withstanding  the  slow 
relaxation  encountered  by  us  and  other  authors  during  (fc  f/  5  T)  measure- 

,  p 

ments,  stress -strain  relaxations  derived  from  force -temperature  data  should 
in  general  be  closer  to  mechanical  equilibrium,  especially  when  the  data  have 
been  derived  after  some  temperature  cycling  at  constant  strain.  This  should 
be  true  for  the  data  on  which  we  based  the  previous  discussion  of  values. 

In  contrast  to  this,  the  "adiabatic"  stress-strain  data  should  be  extensively 
of  non-equilibrium  nature.  The  rapidly  relaxing  stress  peak  must  then  be  due 
to  the  presence  of  relaxation  mechanisms  which  are  slower  than  our  fast¬ 
stretching  experiments.  While  the  data  are  too  few  to  allow  any  other  con¬ 
clusions,  the  principle,  if  established,  should  offer  one  important  approach 
to  the  so  elusive  problem  of  chain  entanglement  .  In  a  slighty  different  way, 
the  force  tracings  during  retraction  can  throw  light  on  the  mechanism  of 
internal  viscosity  and,  thermal  extension  and  the  re-partition  of  energetic 
degrees  of  freedom. 

Finally  some  suggestions  for  further  work  in  the  investigation  of  Rubber - 
like  elasticity  may  be  permissible.  The  nature  of  of  the  Mooney-Rivlin 
equation  has  been  puzzling  for  a  long  time.  It  would  be  valuable  to  study  the 
internal  energy  contribution  to  C  and  to  the  total  retractive  force  by  employing 
systematically  a  number  of  well  defined  elastomers  of  various  different  types. 
To  learn  more  of  the  effect  of  conformational  changes  as  related  to  bond  angles 
and  rotational  energies  upon  the  »ign  and  magnitude  of  fe/f  would  be  of  major 
interest.  The  same  is  true  for  extensive  work  in  dilatometry  and  further 
studies  of  adiabatic  stretching.  The  experimental  system  developed  here 
would  appear  to  have  many  advantages  over  others.  For  adiabatic  work, 
since  no  direct  measurements  of  specific  heat  as  a  function  of  strain  have 
been  reported,  and  it  would  be  most  worthwhile  to  compare  the  values  of  Cp  j_, 
obtained  from  the  measurement  of  force  temperature  and  elongation  curves 
with  those  values  obtained  directly  from  calorimetry. 
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Force-temperature  experiments  at  constant  length,  corrected 
for  thermal  expansion,  allow  to  re-evaluate  and  re-lnterpret  the 
thermoelastic  Inversion  point.  The  extent  of  the  temperature 
dependence  of  the  Mooney-RI vl In  constants  C,  and  C,  could  be  derived 
and  their  values  be  discussed  In  te.ms  of  the  Internal  energy 
contribution. 

Stress-strain  data  from  fast  (near  adiabatic)  rubber  stretch¬ 
ing  show  the  process  to  be  non-equl 1 Ibrlum  In  nature,  and  rate 
dependent  force  peaks  develop.  The  calculated  values  of  C  *  thus 
are  found  to  differ  from  theoretical  expectations.  p* 

\  a  new  dllatometer  was  constructed  which  allows  to  follow 
volume  changes  with  stretching.  For  the  first  time  complications 
due  to  swelling  could  be  overcome.  Extremely  small  volume  Increases 
were  determined  during  the  early  exploratory  runs,  but  the  dllatometer 
has  to  be  reconstructed  to  eliminate  metal  corrosion. 
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